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RESUMO 


O presente trabalho visa o estudo e analise de desempenho de diferentes 
criterios de falha dispomveis na literatura, a diferentes estados triaxiais de tensao. 
Para isso, utiliza-se um corpo de prova ad hoc, aqui denominado Bifailure, que foi 
projetado com a finalidade de se obter falha em estado de tensao proximo ao 
cisalhamento puro e em estado de tensao de alta triaxialidade, em um unico teste de 
tragao. Foram realizados ensaios experimentais e simulagoes numericas com o 
programa comercial de elementos finitos LS-Dyna®. A fim de utilizar criterios de falha 
recentes, ainda nao dispomveis, uma subrotina de elementos finitos foi desenvolvida 
e implementada. Apos as simulagoes, concluiu-se que, dentre os criterios 
analisados, nao houve um criterio que representasse perfeitamente e 
simultaneamente os dois tipos de falha do especime Bifailure. Por fim, o estudo 
concluiu que o grau de acuracidade de um criterio de falha esta relacionado ao 
numero de parametros necessarios e aos testes experimentais que foram realizados 
para sua caracterizagao. 


Palavras-chave: Imperfeigoes e falhas dos materiais. Metodo dos elementos finitos. 



ABSTRACT 


The present work aims the study and analysis of performance of different 
failure criteria available in literature, in different triaxial stress states. For that, it is 
used an ad hoc specimen, here named Bifailure, that was designed with the purpose 
of obtain failure in a stress state near pure shear and in a stress state of high 
triaxiality, in a single tension test. Experimental tests and numerical simulations were 
done with commercial finite element software LS-Dyna®. In order to use recent failure 
criteria, that are still not available, a finite element subroutine was developed and 
implemented. After simulations, it was concluded that, among analysed criteria, there 
was not a criterion that represented perfectly and simultaneously the two types of 
failure of Bifailure specimen. Finally, the study concluded that the degree of accuracy 
of a failure criterion is related to the number of necessary parameters and to the 
experimental tests that were done for its characterization. 

Keywords: Imperfections and failure of materials. Finite element method. 
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1 INTRODUQAO 


A falha corresponde a perda da capacidade de suportar carga por parte de 
uma estrutura ou componente, e e resultado da acumulagao de microdefeitos no 
material, o que constitui o dano. A falha nao resulta necessariamente no colapso 
global de uma estrutura, podendo tambem ser considerada como urn fenomeno 
localizado. 

O fenomeno de falha esta presente no cotidiano das pessoas, em varias 
situagoes onde sao verificadas rupturas em materiais, por exemplo: urn copo de 
vidro quebrando-se ao cair no chao, uma folha de papel sendo rasgada, uma trinca 
em uma estrutura metalica ou de concreto, etc. Sao fenomenos de falha do dia-a- 
dia, alguns controlados e outros nao. 

A falha vem sendo investigada desde os tempos antigos, desde a era das 
cavernas, quando os ancestrais dos seres humanos utilizavam a falha de pedras 
para fabricar utensflios, entre eles armas. 

A revolugao industrial, por sua vez, desencadeou uma longa pesquisa sobre a 
resistencia dos materiais e desenvolvimento de materiais novos. Os metais em 
particular, passaram a integrar estruturas de engenharia como pontes, vefculos 
terrestres, navios, aeronaves, ediffcios, entre outros dispositivos. Assim, o estudo 
sobre a resistencia e ductilidade de metais passou a ser de grande valor na 
engenharia de estruturas. 

Para haver falha ductil e necessaria a ocorrencia de grande quantidade de 
plasticidade, de modo que no cotidiano das pessoas geralmente isso esta presente 
em eventos de impacto como em colisoes de navios com obstaculos, impacto de 
passaros com aeronaves, perfuragoes de blindagens com projeteis e crash-tests de 
vefculos. 

Dessa forma, o estudo de falha e importante para a solugao de problemas 
que ocorrem em varias areas da engenharia. A economia de material em estruturas 
e urn exemplo. Na industria automotiva, na decada de 80, havia urn aumento de 20 
quilos por ano na massa do vefculo e na ultima decada este aumento alcangou 
ainda uma taxa de 10 quilos por ano. 
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Segundo Christlein e Hambrecht (2008), isto ocorre principalmente devido ao 
aumento dos itens de seguranga como o airbag e o antilock braking system (ABS), 
sendo este responsavel por 30% desse aumento de massa. Meihorias no conforto e 
acabamento interno sao responsaveis por 22% e 15%, respectivamente. 

Deste modo, devem-se estudar alternativas para reforgar a estrutura do 
vefculo sem comprometer sua massa. Isso leva a industria em questao a utilizar 
materials cada vez mais leves (Picket et al., 2004) - como compositos, alurmnio, 
magnesio e agos de alta resistencia - que apresentam ductilidade menor que os 
agos comuns, apesar do alto limite de resistencia (Yamashita et al, 2003). Assim, a 
menor ductilidade desses materiais contribui para aumentar a possibilidade de falha. 





Figura 1. Falha em casco de navio apos colisao. Extrafdo do site 

http://www.brighthubengineering.com/seafaring/36193-what-is-ship-collision/ em 22/12/2012. 


O fenomeno de falha tambem ocorre em acidentes com navios (Urban, 2003), 
quando geralmente ocorre ruptura do casco decorrente de colisao. Anualmente 
nesses acidentes perdem-se cerca de 230 embarcagoes, a maioria de pequeno 
porte, e mais de 1000 vidas humanas, aiem de perdas de vida selvagem, quando 
ocorrem acidentes com danos ao meio ambiente (Tornqvist, 2003). Entre os 
acidentes navais onde a falha esta presente podem-se citar, por exemplo, o 
derramamento de oleo, resultante do encalhamento do petroleiro Exxon Valdez, nas 
aguas do Alaska em 1989; o naufragio do navio cruzeiro Estonia em 1994, este 
causado pela falha resultante do estresse gerado quando feita navegagao a alta 
velocidade em mar agitado; e o vazamento de combustfvel do navio maltes Sea Bird 
em 2008, resultante de uma rachadura no tanque de combustfvel apos colisao com 
o navio grego Syrus. 
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A falha e um fenomeno importante a ser considerado tambem no projeto de 
estruturas aeronauticas. A estrutura de um aviao deve ser projetada para resistir ou 
faihar de forma controlada em eventos de coiisoes com passaros e objetos ou com 
fragmentos de turbina em um evento de blade-off, conforme ilustrado na Figura 2. 
Segundo Walvekar (2010), estima-se que 30000 coiisoes de aeronaves com 
passaros ocorrem anualmente, resultando em perdas de 3 bilhoes de dolares para a 
aviagao mundial. 



Figura 2. a) Exemplo de falha em fuselagem de aviao, resultante de colisao com passaro. b) Falha no 
interior de carcaga de motor de aviao, em evento de blade-off. Figuras extratdas, respectivamente, 
dos sites http://gawker.com/5930769/bird-strike-does-terrifying-damage-to-nose-of-united- 
airlines-plane em 22/12/2012 e http://www.nasa.gov/centers/glenn/news/AF/2007/ 
Jan07_printall_prt.htm em 24/12/2012. 


Como aplicagoes militares para o estudo de falha podem ser citadas o projeto 
de blindagens e estruturas que resistam a explosoes. 

Ensaios experimentais de prototipos sao demorados e custosos. Dessa 
forma, a competitividade industrial pelo desenvolvimento de produtos mais eficientes 
leva a crescente necessidade em se reduzir custos e tempo de desenvolvimento de 
produtos. Nesse contexto, a simulagao do comportamento estrutural dos materiais 
pelo Metodo dos Elementos Finitos (MEF) tornou-se uma ferramenta de grande 
valor. Com a constante evolugao dos codigos, assim como da capacidade 
computacional, e possfvel analisar o comportamento da estrutura desde sua fase 
elastica ate a falha em situagoes complexas de geometria e carregamento. 
Exemplos de simulagao numerica de falha atraves do MEF sao ilustrados na Figura 
3. 
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Figura 3. Aplicagoes da simulagao de falha com MEF. a) Impacto em para-choque. Extrafdo de 
Bugelli e Driemeier (2010); b) Impacto em casco de navio. Extrafdo de Tornqvist (2003); c) Impacto de 
passaro em asa de aviao. Extrafdo do site http://www.idacireland.com/autodyn em 22/12/2012. 


A Figura 3a mostra a propagagao de uma falha em simulagao de teste de 
resistencia de para-choque. Na Figura 3b e simulado um evento de colisao de navio, 
com ocorrencia de falha no casco. Por sua vez, a Figura 3c apresenta a simulagao 
de impacto de ave com asa de aviao, na qual verifica-se se a asa falha ou nao. 

Modelos que mimetizem o comportamento real do material, porem, sao 
essenciais para que os resultados numericos reflitam o comportamento real da 
estrutura. Os modelos mais completos disponfveis na literatura estabelecem 
relagoes constitutivas e criterios de falha como fungao de grandes deformagoes, 
efeito de taxas de deformagao e temperatura, perda de rigidez por danejamento, 
evolugao do dano ate surgimento da falha. Porem, ensaios experimentais sao 
sempre muito custosos, e a caracterizagao de alguns desses modelos exige 
equipamentos sofisticados e licengas especiais. 

Atualmente, portanto, projetos devem minimizar o numero de ensaios 
experimentais, de modo a colher com eficiencia dados que caracterizem o modelo 
de material adequado as condigoes de estudo. Esses dados, analisados, fornecem o 
input necessario para simulagao numerica do problema. Dessa forma, podem-se 
racionalizar os investimentos, direcionando novos poucos testes experimentais, em 
situagoes reais, para validagao da resposta numerica. 

O avango da pesquisa relacionada a area possibilita o desenvolvimento de 
modelos constitutivos de material, assim como criterios de falha, que apresentam 
boa acuracidade na representagao dos fenomenos ffsicos modelados. Nota-se 
inclusive que existe uma forte parceria entre empresas e universidades no estudo e 
desenvolvimento de novos modelos de material e de falha (Ockewitz et al, 2006). 
Ainda, ao se analisar a documentagao de programas comerciais de elementos 
finitos, como o LS-Dyna® (2009), verifica-se que necessidades das industrias 
(principalmente a automobilfstica) sao atendidas com relagao a adaptagao e inclusao 
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de novas funcionalidades nos codigos de MEF. Em se tratando do LS-Dyna i®, como 
exemplo, pode-se citar o modelo de previsao de falha em pontos de solda 
desenvolvido pela Toyota e implementado no programa. 

Atuaimente, existem diferentes criterios de falha dispomveis em programas 
comerciais de elementos finitos, como LS-Dyna i®, Abaqus® e Pam-Crash®. No 
entanto, para facilitar a aplicagao pratica, e desejado que o criterio adotado possa 
ser calibrado com urn numero razoavel de testes experimentais simples, e que o 
significado ffsico dos parametros que o caracterizam seja conhecido. 

Diante do exposto, neste trabalho foram analisados diversos modelos de falha 
dispomveis na literatura, atraves do uso do software comercial LS-Dyna®. Esses 
modelos de falha foram implementados no LS-Dyna ® atraves de uma subrotina. O 
modelo de material implementado foi o elastoplastico de von Mises tradicional. Para 
comparagao de acuracidade entre os modelos foi utilizado urn especime ad-hoc 
desenvolvido por Alves e Driemeier (2010), denominado aqui como Bifailure. 

O material em estudo e uma chapa de liga de alurmnio 2024-T351, de 
espessura nominal 9,525mm, similar ao utilizado em Wierzbicki et al. (2005). Foram 
utilizados parametros de material dispomveis na literatura em conjunto com valores 
obtidos de ensaios experimentais de tragao uniaxial. 

Neste trabalho, o Capftulo 2 trata da fundamentagao teorica, com conteudo 
necessario ao entendimento sobre falha como invariantes de tensao e 
microestrutura dos materiais. No Capftulo 3, sao explicados diferentes modelos de 
falha dispomveis na literatura. No Capftulo 4, os testes e resultados experimentais 
sao detalhados. No Capftulo 5 a subrotina de elementos finitos implementada no LS- 
Dyna® e explicada detalhadamente. E no Capftulo 6 os diferentes criterios de falha 
sao aplicados na simulagao dos experimentos com o especime Bifailure, sendo 
discutidos os resultados numericos gerados por cada urn. 
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2 FUNDAMENTOS TEORICOS DO ESTUDO DE FALHA 


A falha estrutural refere-se a perda da capacidade de suportar carga de uma 
estrutura ou componente. Esta se inicia quando o material atinge seu limite de 
resistencia, causando sua separagao em duas ou mais partes. A falha em um 
material e comumente classificada como sendo fragil ou ductil, podendo ocorrer um 
tipo ou outro de falha, ou ambos, de acordo com o tipo de material, o estado de 
tensao, historico de carregamento e temperatura. Os modos de falha no momento 
de ruptura podem variar, entretanto, para materiais ducteis como ago e alummio, o 
fenomeno pode ser ilustrado pela curva tensao-deformagao de um ensaio de tragao 
esquematizada na Figura 4. 

Percebe-se pela figura que em testes de tragao, usados para obter dados 
relativos ao comportamento mecanico do material, a falha e evidente. Distinguem-se 
varias etapas do comportamento mecanico do material, em trechos: 

OA - comportamento linear, com deformagoes nao permanentes; 

AB - plastificagao e encruamento, fenomenos associados ao movimento de 
discordances na rede cristalina; 

BC - danejamento, onde se inicia o rompimento de ligagoes atomicas ate a 
ruptura do material em C. 



O fenomeno de falha em materiais ducteis envolve desde o comportamento 
linear de sua fase elastica ate, depois de atingido o limite elastico do material, sua 
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plastificagao e encruamento (fenomenos associados ao movimento de discordances 
na rede cristalina), seguido de dano (onde ocorre o rompimento de ligagoes 
atomicas). Como consequencia da evolugao do dano, ocorre localizagao de 
deformagoes e consequente estricgao do corpo de prova; o material perde 
resistencia e rigidez; ocorre crescimento e coalescencia de vazios, levando, 
finalmente, ao surgimento de uma macrotrinca que se propaga ate a falha. E, 
portanto, urn fenomeno complicado da engenharia estrutural. 

O corpo de prova de ago inoxidavel (austemtico) da Figura 5 foi testado em 
urn ensaio de tragao quase-estatico. Percebe-se pela figura que ha consideravei 
plastificagao e estricgao em uma pequena regiao do corpo de prova, onde ocorre a 
ruptura. Essa estricgao esta relacionada com o imcio do danejamento do material e 
ocorre no ponto B da Figura 5. Nessa regiao, o estado de tensao nao e mais 
uniaxial, podendo ser medido experimentalmente atraves de tecnicas especiais de 
tratamento de imagens (como o Digital Image Correlation - DIC) ou, simplesmente, 
no caso de corpos de prova cilfndricos, atraves do metodo de Bridgman (como 
descrito em Ling, 1996). 


(b) 

(a) 


Figura 5. Corpo de prova de ago (a) antes e (b) apos ensaio de tragao uniaxial. Teste realizado no 
laboratorio do Grupo de Mecanica dos Solidos e Impacto em Estruturas (GMSIE). 

Varias formas de tratar o fenomeno de falha foram desenvolvidas, entre elas 
modelos de falha abrupta, modelos micromecanicos e mecanica do dano contfnuo, 
que serao explicados neste trabalho. 

O fenomeno de falha geralmente envolve nao linearidades de material e 
geometrica, assim torna-se conveniente utilizar metodos numericos para prever a 
falha, sendo que com o constante aumento da capacidade de processamento dos 
computadores, os codigos de elementos finitos podem se tornar capazes de lidar 
com modelos mais complexos de falha. Por outro lado, a praticidade de calibragao 
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de um criterio de falha e uma prioridade na definigao de um criterio de falha robusto, 
que tambem precisa ser capaz de ser implementado numericamente e prever com 
boa precisao a falha em diferentes situagoes a fim de poder ser aplicado com 
facilidade na industria. 

Para o metodo numerico ser utilizado corretamente, e necessario conhecer 
bem o fenomeno estudado, que tern suas explicagoes na microestrutura. Alem disso, 
e necessario modelar o material e a falha de forma consistente. 

De acordo com Wu et al. (2010), o comportamento estrutural de um material 
depende das propriedades do material, da microestrutura e condigoes de 
carregamento (tensoes, deformagoes e temperatura). 


2.1 Microestrutura 


Quando um metal solidifica-se, forma-se uma estrutura policristalina na qual 
cada cristal e chamado de grao e cada grao possui uma orientagao cristalografica 
definida. As fronteiras entre os graos sao chamadas de contornos de grao, conforme 
ilustrado na Figura 6. 



Figura 6. Estrutura granular de ago, mostrando nitidamente os contornos de grao. Extraido do site 

http://www.arcelormittal.com/automotive/saturnus/sheets/catalogue.pl?id_sheet=l5&header=&l 

anguage=EN em 30/12/2011. 

As propriedades mecanicas da estrutura policristalina sao uma media das 
propriedades mecanicas de cada grao anisotropico, assim como um todo a estrutura 
comporta-se como sendo isotropica. 
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A estrutura atomica nao e perfeita: nela existem defeitos como lacunas, 
atomos substitucionais ou intersticiais, discordancias ou vazios. A plasticidade e 
embasada no deslocamento de discordancias na rede cristalina, sendo que uma 
discordancia e um defeito linear ou unidimensional que gera uma distorgao local da 
rede cristalina, conforme Figura 7. 



Figura 7. Discordancia do tipo aresta, indicada pela seta, gerando distorgao na rede cristalina. 

Adaptado de Callister (2007). 

O deslocamento de discordancias, induzido por tensoes cisalhantes, ocorre 
ao longo de pianos cristalograficos chamados de pianos de escorregamento. 

Quanto maior a dureza e resistencia de um material, menor e a facilidade com 
que as discordancias conseguem se mover no reticulo cristalino. Durante o processo 
de deformagao plastica o numero de discordancias aumenta consideravelmente, 
geradas pela muitiplicagao de discordancias pre-existentes e por concentradores de 
tensao como contornos de grao e microdefeitos do material. Como a interagao entre 
duas discordancias e de ordem repulsiva, quando as discordancias interagem entre 
si o movimento de uma discordancia acaba sendo dificultado pela presenga de 
outras discordancias, e com isso temos o aumento da resistencia mecanica do 
metal. 

Em um ensaio de tragao, ao se atingir o limite de resistencia a tragao do 
material, a deformagao plastica torna-se nao homogenea, concentrando-se em uma 
pequena porgao do corpo de prova, e gerando uma zona de estrangulamento da 
segao conhecida por estricgao. A deformagao torna-se nao homogenea pelo fato de 
a diminuigao da segao transversal ocorrer mais rapidamente que o encruamento. A 
seguir, ocorre a ruptura (falha) do material. 

Em falha ductil, que e o foco deste trabalho, ha consideravel deformagao 
plastica antes da ruptura. Geralmente quanto maior a pureza de um material ductil, 
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maior a deformagao piastica suportada ate a ruptura. O termo “falha ductil” 
geralmente esta relacionado com a falha por nucleagao, crescimento e propagagao 
de vazios (Ballard, 1997). 

A influencia de temperatura e evidente na definigao entre ocorrencia de falha 
mais ductil ou mais fragil, apesar do comportamento ductil do material a temperatura 
ambiente. Quanto maior a temperatura, mais intenso e o movimento de vibragao dos 
atomos. Assim maior e a facilidade de restabelecimento de ligagoes atomicas 
quando o material sofre plastificagao, pela redugao da resistencia ao movimento 
relativo entre atomos. Dessa forma, a falha no material ocorre com presenga de 
grandes deformagoes plasticas, o que e caracteristico da falha ductil. Quando a 
temperatura e reduzida ocorre o efeito inverso: o movimento de atomos ao longo da 
rede cristalina e dificultado, assim como o restabelecimento de ligagoes atomicas 
rompidas. Quando a falha ocorre, a dificuldade de restabelecimento de ligagoes 
atomicas dificulta a progressao da deformagao piastica, o que caracteriza o modo de 
falha fragil. 


2.2 Invariantes de tensao 


Considerando-se as componentes do tensor de tensoes a ij e as tensoes 
principais a x , <j 2 e define-se a tensao hidrostatica como: 


1 . 1 If \ 

= 3 I] = 2 a " = 2^ +(7 i +(7 i’ 


( 1 ) 


onde /, representa o primeiro invariante do tensor de tensoes e a notagao com 
indices iguais representa urn somatorio. Considerando Sjj como as componentes do 
tensor desviador, e Si, S2e S3 os seus valores principais, tem-se: 


(2) 
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onde 8 tj e o delta de Kronecker. Em forma matricial, 8 tj corresponderia a matriz 
identidade. O segundo e o terceiro invariantes do tensor desviador sao: 



(3) 


(4) 


A tensao equivalente de von Mises (<x) pode ser expressa por: 


g = sj3J 2 


(5) 


Considere urn volume diferencial submetido ao estado de tensoes 
tridimensional no espago das tensoes principals (g 1 ,g 2 ,g 3 ). Na Figura 8a, a reta 

paralela a ON, ou reta h , que passa pela origem com angulos iguais entre os tres 
eixos e chamada de eixo hidrostatico, e cada ponto desse eixo corresponde a urn 
estado hidrostatico de tensoes. O piano perpendicular a reta h que passa pela 
origem e chamado de piano n, piano octaedrico, ou ainda piano desviador, e nesse 
piano a tensao hidrostatica e nuia. Considerando urn ponto P no espago de tensoes 
com componentes {g v g 2 ,g 3 ), o vetor OP pode ser decomposto em uma componente 

ON, com norma a, paralela ao eixo h , relacionada a tensao hidrostatica, e em uma 

componente NP, com norma r, perpendicular a h , relacionada a parte desviadora, tal 
que: 


a = 


V3 a, 


( 6 ) 



(7) 


Desta forma, o fator de triaxiaiidade e definido como sendo: 


5 !L= V2 a 

G 3 r 


(8) 
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Figura 8. (a) Cone representative de um unico estado de triaxialidade, (b) Projegao de um estado de 

tensoes no piano n 

A co-tangente do angulo </> da Figura 8a e proporcional ao fator de 
triaxialidade. Portanto demonstra-se que para um determinado fator de triaxialidade 
existe um infinito numero de estados de tensao, representados pela superffeie do 
“cone” da Figura 8a. Ainda, para completar a determinagao do estado de tensao e 
necessaria a inclusao do angulo de Lode, 0 L , o qual numericamente representa a 
magnitude da tensao principal intermediaria <7, em relagao a maior e a menor 
componente de tensao <7, e (t 3 , e geometricamente relaciona-se com a posigao do 

ponto P na circunferencia representada na Figura 8b. 

As tensoes desviadoras podem ser obtidas pela projegao do vetor de tensoes 
OP no piano k , conforme ilustrado na Figura 8b. Por essa razao que o piano n 
tambem e chamado de piano desviador. 


o 2 



Figura 9. Projegao do vetor de tensoes no piano n. Adaptado de Xue (2007). 
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O angulo de Lode corresponde ao angulo 0 L da Figura 9. Na Figura 9, alem 
do angulo de Lode 6 ,, aparecem as relagoes entre as tensoes principals no piano 

n. 

O piano k aparece esquematizado na Figura 9. Nesse piano, o angulo entre 
os eixos das tensoes principals e igual a 120 Q . Desse modo, tem-se que: 


tan 0 L 


(7 2 -(J x sen30° - <7 3 sen30° 
<J , cos 30° - <J 3 cos 30° 


0) 


Simplificando, obtem-se a definigao do angulo de Lode 0 L : 


0, 


= tan 


[ 1 2<j 2 — (Tj — (X 3 

V3 cr.-cr. 


( 10 ) 


tal que — <0 L < — . Dessa forma, tem-se que o angulo de Lode, para alguns 

6 6 

estados de tensao mais conhecidos, e definido como: 

a x ± 0, <x, = <r 3 = 0 com 0 L = tensao uniaxial 

<y, = —<t 3 , o '2 — 0 com d L =0 cisalhamento 

<7 l =<7 2 , <r 3 <0 com d L = K /§ compressao uniaxial 

A relagao entre 0 L e o terceiro invariante 7 3 e dada por: 

f=^-4=“> s W on 

2 (X 

onde £ e o terceiro invariante da parte desviadora normalizado, de forma que 
-1<£<1. O parametro E, caracteriza a posigao (magnitude) da segunda tensao 
principal o 2 em relagao a primeira e terceira tensoes principals a x e a 3 . Quando 
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£ = - 1 , tem-se <j x = <r 2 > <j 3 ; quando £ = 1 , tem-se <t, > a 2 = o 3 ; enquanto que £ = 0 
para <r 2 = (<x, + a 3 )/2. 

Quando o angulo de Lode e medido no sentido anti-horario, a partir da linha 
de ( 7 , na Figura 9, variando assim de zero ate tu/ 3, utiliza-se tambem a 
normalizagao do angulo de Lode, conhecida como parametro de Lode (Bai e 
Wierzbicki, 2008): 


a , 6(9, 2 - 

<9 =1-- = 1- acosg 

71 K 


( 12 ) 


tal que — 1 < 6 > < 1 . Considerando-se o parametro de Lode 6 , tem-se que <9=1 
corresponde a tragao axissimetrica, 0=0 corresponde a cisalhamento, e <9=-l 
corresponde a compressao axissimetrica ou estado de tragao equi-biaxial. 

O estado piano de tensao {cr 3 = 0 ) univocamente relaciona 77 com £ ou 0 : 


£ = cos( 3# l ) = cos 




(1-0) 
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: -7 

2 1 


^ 2 o 

7 — 

V 3y 


(13) 


A equagao anterior apresenta tres raizes, uma correspondendo a 
cisalhamento puro (77 = 0 , 7 = 0 ) e as outras duas ao estado piano de deformagao 

( t 7 = ±i/V3,7 = 0 ). 

Em smtese, urn estado de tensao representado pelas tensoes principals no 
sistema cartesiano de coordenadas (<7 1 ,<j 2 ,<j 3 ), tambem pode ser representado no 

sistema cilindrico ( o m ,o,0 L ) ou no sistema esferico de coordenadas (a,jj,0 L ). 

Outro parametro derivado de invariantes e que e utilizado em trabalhos de 
falha e 0 produto de triaxialidade: 


* 



_ _ J3 


(14) 


A maioria dos criterios de falha utiliza os invariantes e parametros citados 
acima na definigao de fungoes de deformagao de falha e c . 
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3 REVISAO DA LITERATURA 


O fenomeno de faiha nos metais e resultado de um processo complexo que 
ocorre em escala micrometrica. Tal processo envolve a evolugao e acumulagao de 
defeitos microscopicos como microvazios e bandas de cisalhamento. Esses defeitos 
constituem o chamado “dano” do material, afetando as propriedades mecanicas do 
mesmo. 

No passado, varias abordagens foram propostas para o tratamento de faiha 
ductil, como a mecanica da fratura, modelo de crescimento de vazios, modelo de 
poroplasticidade, mecanica do dano e varios modelos empmcos. A mecanica de 
fratura lida com a propagagao de uma trinca, o que e diferente do tema tratado nesta 
dissertagao. Assim, por esse motivo, a trinca originada pelo crescimento e 
propagagao de microvazios no material sera chamada aqui de faiha e nao de fratura. 

As observagoes experimentais sobre faiha ductil em geral mostram que (Xue 
e Wierzbicki, 2008): 1) a aplicagao de pressao aumenta a ductilidade do material; 2) 
a ductilidade do material em estado piano de deformagao (ou cisalhamento) e menor 
do que a ductilidade do mesmo sob o estado de tragao ou compressao quando a 
pressao mantem-se constante; 3) o dano no material se acumula de forma 
acelerada, inicialmente a uma baixa taxa e posteriormente a uma taxa mais alta a 
medida que a faiha se aproxima; 4) o dano do material afeta a sua resistencia 
mecanica. 

A revisao da literatura cientffica dispomvel revela grande empenho dos 
centros de pesquisa no estudo de criterios de faiha mais representatives dos 
fenomenos analisados, considerando a viabilidade de calibragao desses criterios 
atraves da minimizagao de testes experimentais. Apesar disso, deve-se ressaltar 
que a aplicagao inadequada desses modelos para diversos estados complexos de 
deformagao pode resultar em previsoes erroneas. 

Inicialmente, surgiram modelos de faiha que consideravam a acumulagao de 
microdefeitos no material. Considerou-se que o surgimento, crescimento e 
propagagao de micro-vazios era responsavel pela faiha ductil dos materiais. A 
primeira analise microscopica de faiha foi proposta por McClintock et al. (1966) que 
assumiu uma distribuigao de vazios cilmdricos no material e estudou a evolugao do 
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crescimento desses vazios. Posteriormente, Rice e Tracey (1969) estudaram o 
crescimento e mudanga de forma de vazios esfericos nos materiais quando sujeitos 
a carregamento uniaxial. Os autores propuseram que o aumento de volume era mais 
significativo do que a mudanga de forma dos vazios e portanto escreveram a 
deformagao de falha como fungao do fator de triaxialidade, que e a razao entre a 
tensao hidrostatica e tensao equivalente. No entanto, nos dois modelos anteriores 
nao foi estudada a interagao entre microvazios no material. Assim, Leroy et al. 
(1981) modificou o modelo de Rice e Tracey para incluir o efeito de nucleagao e 
interagao entre vazios. Todavia, esses modelos nao consideram a plasticidade como 
acoplada ao dano acumulado pelo material. De acordo com esses modelos, a falha 
ocorre quando uma variavel interna de dano atinge urn valor crftico, resultando em 
uma perda repentina da capacidade de suportar carga do material. Dessa forma, 
surgiram modelos com acoplamento entre plasticidade e dano, como o proposto por 
Gurson (1977), que considera o material como sendo urn meio poroso, no qual urn 
fator denominado fragao volumetrica de vazios controla a ocorrencia de falha. Aqui a 
falha e resultante da coalescencia de vazios, e ocorre quando a fragao volumetrica 
de vazios se iguala a unidade. O modelo de Gurson nao e realfstico, pois, por 
exemplo, nao e capaz de representar a falha sob condigoes de cisalhamento, onde o 
mecanismo de crescimento de vazios e inativo. Desse modo, o modelo original de 
Gurson foi repetidamente modificado para levar em consideragao processos 
adicionais responsaveis pela deterioragao do material e consequente falha ductil. 
Por exemplo, a fim de aprimorar o modelo de Gurson, Tvergaard (1981) e Tvergaard 
e Needleman (1984) propoem considerar a fragao volumetrica de vazios uma fungao 
da nucleagao de novos vazios e do crescimento de vazios existentes. No entanto, 
apesar dos aprimoramentos do modelo de Gurson, o grande numero de parametros 
do modelo e o fato destes serem fortemente acoplados tornou diffcil a calibragao e 
aplicagao do modelo. Outro tipo de modelo de falha acoplado, denominado de 
Mecanica do Dano Contfnuo, foi proposto por Lemaitre (1985). Aqui a degradagao 
do material e modelada a partir de uma variavel interna de dano enquanto que as 
equagoes constitutivas sao derivadas da primeira e segunda leis da termodinamica 
para meios contfnuos. 

Para aplicagoes industriais, sao utilizados modelos de falha fenomenologicos 
empfricos. Apesar de nao possufrem embasamento micromecanico, a simpiicidade e 
acuracidade destes modelos os tornam mais aceitos em aplicagoes de engenharia. 
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O criterio de falha empirico mais simples e o baseado na maxima deformagao 
plastica equivalente. No entanto, demonstra-se que a deformagao equivalente por si 
so nao e suficiente para indicar a ocorrencia de falha (Alves e Jones, 1999). 

Atualmente, diversos autores confirmam que a falha e fortemente dependente 
do estado de triaxialidade do material, alem do mvel de deformagao plastica [Alves e 
Jones (1999), Bao e Wierzbicki (2004), Borvik et al (2004), Wierzbicki et al (2005), 
Bao e Wierzbicki (2005), Ockewitz et al (2006) e Grytten et al (2009a)]. O modelo de 
falha de Johnson-Cook (1985), que tambem e fungao da triaxialidade, e bastante 
conhecido por ser o primeiro modelo de falha que tambem e fungao da taxa de 
deformagao e temperatura. 

Para a caiibragao dos criterios baseados no fator de triaxialidade e necessario 
gerar dados experimentais de tensao, triaxialidade e deformagao equivalente no 
instante de falha. Assim, a forma mais comum de se obter diferentes faixas de 
triaxialidade em experimentos era conduzir testes de tragao em especimes 
ciiindricos sem e com diferentes entalhes. A seguir, os tensores de tensao e 
deformagao no instante de falha desses especimes poderiam ser estimados por 
metodos como o de Bridgman (1956). Porem, metodos como o de Bridgman apenas 
poderiam gerar dados para altas triaxialidades, limitando assim a caiibragao e 
aplicagao dos modelos de falha para altas triaxialidades. Desse modo, a posterior 
introdugao das analises pelo Metodo dos Elementos Finitos foi essencial para prever 
os valores de tensao e deformagao na regiao de falha, permitindo uma melhor 
caiibragao dos modelos de falha. 

Kim et al. (2003, 2004, 2007) e Bai e Wierzbicki (2008) demonstraram que 
apenas o fator de triaxialidade nao e suficiente para caracterizar o estado de tensoes 
de falha, sendo necessario considerar tambem urn novo parametro, o angulo de 
Lode. Este fato e confirmado pela observagao experimental de que em comparagao 
com os testes de tensao e compressao uniaxiais, uma menor ductilidade e 
observada no teste de torgao (Johnson e Cook, 1985). Desse modo, o criterio de 
Xue-Wierzbicki proposto em Wierzbicki et al. (2005) foi criado para incorporar esta 
observagao experimental, possuindo uma superffcie de deformagao de falha que e 
fungao do fator de triaxialidade e parametro de Lode. Alem disso, investigagoes 
numericas de urn vazio esferico contido em uma celula cubica submetida a tipos 
diferentes de carregamento triaxial (Gao e Kim, 2006) mostraram que o angulo de 
Lode tern uma forte influencia na evolugao do formato do vazio e no processo de 
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coalescencia de vazios. E a influencia do angulo de Lode na ductilidade de metais 
tambem foi demonstrada experimentalmente por Barsoum e Faleskog (2007), 
atraves de ensaios experimentais de um especime cilfndrico especial sob 
carregamentos combinados de tragao e torgao, que produziam valores diferentes de 
triaxialidade e parametro de Lode. 

Os trabalhos cientfficos que estudaram a influencia do angulo de Lode sao 
relativamente recentes, de forma que foi determinada a influencia do angulo de Lode 
tanto em falha quanto na superffcie de escoamento, para diferentes modelos de 
falha e plasticidade disponfveis na literatura. Por exemplo, em Bardet (1990) foi 
estudada a influencia do angulo de Lode no formato da superffcie de escoamento, e 
determinado que o modelo de plasticidade de Drucker-Prager assim como os 
modelos de Tresca e Mohr-Coulomb tern dependencia do angulo de Lode. Alem 
disso, alguns autores sugeriram a inclusao da dependencia do angulo de Lode na 
formulagao de modelos constitutivos tradicionais baseados na superffcie de 
escoamento de von Mises assim como em modelos de evolugao de dano. 
Especificamente, Brunig et al. (2008) e Bai e Wierzbicki (2008) propuseram modelos 
de plasticidade que incluem a influencia dos tres invariantes de tensao na superffcie 
de escoamento. 

Os modelos de plasticidade baseados na tensao de Von Mises ou criterio de 
Tresca podem prever satisfatoriamente o infcio do escoamento de materiais ducteis, 
no entanto e o angulo de Lode que permite representar falha em pianos inclinados, 
como ocorre em testes de tragao ou compressao uniaxial, ou diferengas entre as 
curvas de carregamento de testes de tragao e cisalhamento de um material 
(Chocron et al., 2011). Isto decorre pelo fato de que quando a superffcie de 
escoamento e escrita em fungao do angulo de Lode, a mesma passa a possuir um 
formato nao convexo, adotando a forma de uma “estrela de seis pontas”. 


3.1 Lugar geometrico de falha 


Estudos experimentais sobre falha em diferentes estados de tensoes foram 
conduzidos de forma a se determinar os fatores que influenciam o fenomeno de 
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faiha, de modo que para os modelos empfricos verificou-se ao longo do tempo que a 
falha era dependente da tensao hidrostatica e tensao de von Mises. Esses dois 
fatores, que estao ligados aos invariantes de tensao, foram combinados no fator de 
triaxialidade. A tensao hidrostatica positiva acelera o processo de nucleagao- 
crescimento-coalescencia de vazios e o movimento de bandas de cisalhamento, 
enquanto que a tensao hidrostatica de natureza compressiva tem por efeito mitigar 
esses processos. 

A seguir, verificou-se que a faiha nao era dependente apenas do fator de 
triaxialidade, mas tambem do angulo de Lode, o quai esta relacionado com o terceiro 
invariante de tensao do tensor desviador. No entanto, a sensibilidade da deformagao 
de faiha em relagao ao fator de triaxialidade e ao angulo de Lode e diferente para 
diferentes materiais. Por exemplo, na Figura 10 aparece esquematizada a superffcie 
de deformagao de faiha para dois materiais, o alummio 2024-T351 e o ago 1045, em 
fungao do fator de triaxialidade e angulo de Lode. Percebe-se que a faiha do 
alummio 2024-T351 e fortemente influenciada tanto pelo fator de triaxialidade quanto 
pelo angulo de Lode, enquanto que para o ago 1045 a superffcie de faiha tem baixa 
influencia do angulo de Lode. 



Stress Triaxiality 


Lode angle 


Figura 10. Deformagao de faiha em fungao do 
alumfnio 2024-T351; b) ago 1045 . 


Stress Triaxiality 


Lode angle 


fator de triaxialidade e angulo de Lode para a) 
Extrafdo de Malcher et al. (2012). 


Os mecanismos principais de degradagao de material consistem no 
crescimento de vazios esfericos (influenciado por efeitos de tensao hidrostatica) e 
alongamento de vazios (gerado por efeitos de cisalhamento). 

A Figura 11 esquematiza os parametros de triaxialidade e angulo de Lode que 
podem ser obtidos atraves de especimes comumente utilizados em ensaios 
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experimentais. Na Figura 11, a area “A” representa a regiao na qual a maior parte da 
degradagao do material ocorre por efeitos de cisalhamento, onde ambos os valores 
de triaxialidade e angulo de Lode estao ao redor de zero, assim o mecanismo 
principal de degradagao consiste no alongamento de vazios. Na area “B”, ainda 
existe forte influencia dos efeitos de cisalhamento, observada em condigoes de 
compressao-cisalhamento e tragao-cisalhamento, onde tanto ocorrem mecanismos 
de crescimento de vazios esfericos quanto de alongamento de vazios. Na area “C”, 
os efeitos de cisalhamento sao desprezfveis e o mecanismo de degradagao 
predominante consiste no crescimento de vazios esfericos. 



Figura 11. Representagao de estados de tensao no espago de triaxialidade versus parametro de 

Lode. Adaptado de Malcher et al. (2012). 


Na Figura 11, os valores dos eixos correspondem a valores medios do fator 
de triaxialidade e parametro de Lode, dados por: 


V m = — \v(£)d£ p (15) 

£ f 0 

e. =y-j S(e)de r 

£ f 0 


( 16 ) 
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onde e f corresponde a deformagao plastica equivalente de falha. 

A partir de observagoes experimentais, verifica-se que a falha ductil 
compreende tres fases: acumulo de dano, imcio da falha e propagagao de trinca. O 
infcio da falha ocorre a partir de acumulo de dano, de forma que microscopicamente, 
tal dano e associado com a nucleagao, crescimento, e coalescencia de microvazios, 
movimento de bandas de cisalhamento e propagagao de microtrincas. 
Macroscopicamente, o dano no material e verificado pela redugao da resistencia 
mecanica e ductilidade do mesmo. Essas mudangas no material sao utilizadas para 
a previsao de falha, atraves do controle da verificagao de valores indicativos 
instantaneos ou valores indicativos de falha que variam de forma cumulativa com o 
tempo. Em mecanica do contfnuo, uma variavel indicativa de danejamento no 
material comumente utilizada chama-se “dano”. Apesar de o dano possuir natureza 
anisotropica, por praticidade o mesmo e modelado como sendo isotropico e como 
uma grandeza escalar, o que tern produzido bons resultados. Como o dano e uma 
variavel interna do material que nao pode ser medida diretamente, e necessario criar 
uma correlagao entre o dano e outras variaveis que podem ser mensuraveis. Deste 
modo, os modelos de dano criam a correlagao necessaria para se poder quantificar 
o dano. 

Dois enfoques principals sao dados para a modelagem do comportamento 
mecanico dos materials: o enfoque macroscopico, em que o material e considerado 
em forma agregada, e o enfoque microscopico, em que a microestrutura do material 
e levada em consideragao. 

Modelos macroscopicos sao construidos com base em variaveis de estado 
mensuraveis, como tensores de tensao e deformagao, e as vezes temperatura e 
taxa de deformagao. Demonstra-se experimentalmente que o historico de 
carregamento influencia diretamente a ocorrencia de falha. 

Varios modelos de falha consideram que o dano nao possui influencia na 
plasticidade do material, ou seja, sao modelos desacoplados. Uma vantagem dessa 
abordagem e que a evolugao do dano e a evolugao da plasticidade sao 
desacopladas ao longo do calculo das tensoes no metodo de integragao numerico, 
reduzindo o processamento computacional necessario. Consequentemente, no 
processo de verificagao de falha, a variavel de dano torna-se uma variavel de pos- 
processamento. Por outro lado, o fato de o modelo de dano estar desacoplado 
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resulta em uma queda brusca da curva de carregamento do material quando da 
ocorrencia de falha, resultando em uma representagao nao realfstica dos modos de 
falha, incluindo a falta da perda de resistencia do material proximo do momento de 
falha com a representagao de urn aumento de resistencia nos momentos finais de 
plasticidade que na verdade nao existe. 

Em mecanica do dano contfnuo, considera-se que o dano possui influencia na 
resistencia mecanica do material desde o infcio do carregamento ate a ruptura do 
material. Esta forma de modelagem por outro lado exige maior capacidade de 
processamento computacional. 

Os modelos de dano cumulativos consideram que a falha ocorre quando a 
variavel de dano excede urn valor crftico D c . Normalmente a variavel de dano e 

normalizada em relagao ao valor crftico D c , de forma que o valor crftico de dano a 

ser atingido passa a ser unitario. A equagao que define o incremento do dano 
cumulativo possui a forma geral: 


^— = F(e,T,e) (17) 

de ,. 

onde de p e o incremento de deformagao plastica ao longo dos ciclos de integragao 

no tempo. O lado direito da equagao (17) e fungao de variaveis de estado como 
deformagao, temperatura e taxa de deformagao. Neste trabalho, sao estudados 
criterios de falha que nao consideram a influencia de temperatura nem da taxa de 
deformagao no lado direito da equagao anterior, assim F torna-se F = F(e). 


3.2 Mecanismos fisicos para ocorrencia de falha ductil 


Como ja discutido anteriormente, alem da microestrutura do material e do 
historico de tensoes e deformagoes, urn dos principals fatores que influenciam a 
ocorrencia de falha ductil e a triaxialidade. Assim, sao apresentados a seguir 
diferentes formas de falha ductil, de acordo com a faixa de triaxialidade. 
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Seguindo a ideia de Bao e Wierzbicki (2004, 2005), ilustrada na Figura 12, e 
conveniente distinguir tres estados que levam a falhas distintas, sendo eles: falha 
por crescimento de vazios (alta triaxialidade), falha por cisalhamento (baixa ou nula 
triaxialidade), e a transigao entre os dois (triaxialidade media). Bridgman (1956), por 
exemplo, ja havia definido que quando a tensao hidrostatica e compressiva, a 
ocorrencia do fenomeno de nucleagao-crescimento-coalescencia de vazios e 
dificultada, aumentando assim a ductilidade em metais. 


Equivalent strain to fracture f 



Figura 12. Relagao da deformagao equivalente de falha com a triaxialidade para alummio 2024-T351 

(Bao e Wierzbicki, 2004) 


3.2.1 Alta triaxialidade 


Este modo de falha se caracteriza por gerar superficies asperas e rugosas 
(com alveolos). O infcio do processo ocorre com a nucleagao de vazios nas regioes 
do material com inclusoes de maior tamanho (partfculas de elementos de liga), 
grandes defeitos da rede cristalina ou em regioes com bandas de cisalhamento. Urn 
exemplo que ilustra a nucleagao em partfculas de liga e ilustrado na Figura 13. 

Apos a nucleagao, ocorre o crescimento e propagagao do vazio, sendo que a 
forma de crescimento e mudanga de formato de vazios depende do historico de 
tensoes e deformagoes do material. 
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Superficie lis 
particula se 
destacou da 
matriz cristalir 



Particulas que falharam 


Superficie rugosa: 
faiha ductil da 
matriz cristalina 


5 pm 

Figura 13. Nucleagao de vazios em particulas, e faiha ductil da matriz cristalina para liga de alummio 

2024-T351. Extrafdo de Bao e Treitler (2004). 

Alem da nucleagao e crescimento de vazios, existem diferentes mecanismos 
para a propagagao dos mesmos. Os vazios podem propagar-se devido a estricgao 
da rede cristalina entre dois vazios (Figura 14a), quando dois vazios entram em 
contato (Figura 14b), ou pela ligagao de microvazios entre dois vazios maiores 
(Figura 14c). 


Crescimento de vazios 


aJ A I (A Al A A A . ■ 



Figura 14. (a) Ligagao de vazios por empescogamento; (b) Ligagao de vazios por contato; (c) Ligagao 

de vazios por ligagao de microvazios 


3.2.2 Triaxialidade negativa 

Nao existem estudos detalhados sobre faiha em triaxialidades negativas. No 
entanto, verifica-se que em testes de compressao de cilindros (nos quais a faixa de 
triaxialidade esta entre -0,33 e -0,05) a faiha ocorre nos pianos de tensao cisalhante 
maxima (bandas de cisalhamento), a 45 Q da diregao de carregamento, gerando uma 
superficie de faiha que e relativamente lisa. Podem-se verificar essas caracterfsticas 
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na compressao de ciiindros (Bai e Wierzbicki, 2008), de acordo com a Figura 15. No 
trabalho de Bao e Wierzbicki (2005), inclusive, foi proposto um valor de cut-off de 
-1/3 para a triaxialidade, abaixo do qual nao ha falha. 



Figura 15. Falha em teste de compressao de ciiindros, nos pianos de tensao cisalhante maxima. 

Extraido de Bai e Wierzbicki (2008). 

3.2.3 Falha em triaxialidade intermediaria 


French e Weinrich (1975) estudaram a influencia da pressao em testes de 
tragao de especimes de cobre. Foi verificado que sem aplicagao de pressao, a 
superficie de falha mostra que a mesma ocorreu totalmente por nucleagao, 
crescimento e ligagao de vazios. Por outro lado, a aplicagao de pressao implica na 
redugao da rugosidade da superficie de falha, isto e, ocorre transigao gradual para o 
modo de falha verificado em triaxialidades negativas. Assim, pode-se concluir que na 
faixa de triaxialidades intermediaries, o mecanismo de falha e uma combinagao de 
nucleagao, crescimento e ligagao de vazios e falha nas bandas de cisalhamento. 


3.3 Modelos de falha da mecanica do continuo 


Varios modelos de falha da mecanica do continuo consideram o processo de 
danejamento do material como desacoplado do comportamento plastico do material. 
Assim, nesses modelos a variavel de falha tern por finalidade identificar a sua 
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ocorrencia, nao infiuenciando na perda gradual de resistencia do material ate atingir 
a ruptura. Dessa forma, esses modelos representam uma perda subita da 
capacidade do material de transmitir carga no instante de falha. Por outro lado, a 
mecanica do dano contfnuo considera que o danejamento do material afeta sua 
resistencia. Deste modo, o material perde gradativamente sua resistencia ate a 
ocorrencia de falha. Este e o caso do modelo de dano de Lemaitre. 


3.3.1 Criterio de falha da maxima deformagao plastica equivalente 


Este criterio identifica a falha quando e atingida uma deformagao plastica 
equivalente de ruptura £ /. 



Para urn material incompressfvei s e definido como: 


e = 



(18) 


(19) 


onde £ i> £ 2 , £3 sao as deformagoes plasticas principais. Este criterio e valido para 
todos os estados de tensao, de forma que apesar de ter sido urn dos primeiros 
criterios de falha que surgiu, o mesmo e adotado ainda hoje, pela facilidade em 
identificar a falha na fase de pos-processamento (basta aplicar escalas de cores 
delimitadas pelo valor crftico de deformagao) e pela grande disponibilidade de dados 
de deformagao de ruptura que podem ser encontrados na literatura. 
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3.3.2 Criterio de falha da maxima tensao cisalhante 


Existem inumeras evidencias experimentais de que a falha ocorre no piano de 
tensao cisalhante maxima, como mostrado por Bai e Wierzbicki (2008). A condigao 
de falha e definida como: 


T — (t ) 

max V max / f 


^max = maX i 


0\ <7 2 


<J 2 o 2 


<J 3 (J, 


( 20 ) 

( 21 ) 


onde <7,, <j 2 e <x 3 sao as tensoes principais. 


3.3.3 Criterio de falha de Johnson-Cook 


O criterio de falha de Johnson-Cook possui uma fungao de dano cumulativo 
da forma: 


D = \ 



( 22 ) 


A falha e detectada quando a variavel de dano D atinge o valor unitario. A fungao de 
deformagao de falha, s c, e dada por: 


e c = [D, + D 2 cxp(/J f /;)] 


1 + D 4 log — 

£, 


1 + ZX 


r T-T 0 ^ 


(23) 


onde q e o fator de triaxialidade, s o e a taxa de deformagao de referencia, T f e a 
temperatura de fusao do material, e T 0 e a temperatura de referencia. Neste 



46 


trabalho, por praticidade, o criterio de falha de Johnson-Cook utilizado possui a 
formula: 


£ c = d i+ d i exp {D 3 tj) (24) 

Deste modo, os efeitos da taxa de deformagao e temperatura sobre a 
deformagao de falha foram ignorados. Isso porque para o alummio 2024-T351 os 
efeitos de taxa de deformagao e temperatura sobre a falha sao significativos apenas 
na simulagao de impacto a altas velocidades. Conclusoes similares tambem foram 
tiradas para outros metais, como em Borvik et al. (2001), Borvik et al. (2003) e 
Hopperstad et al. (2003). 

O criterio de falha de Johnson-Cook e amplamente utilizado devido ao fato de 
ser pratico de calibrar e de existir na literatura urn vasto banco de dados sobre os 
parametros de falha para diversos materials (Johnson e Holmquist, 1989). 
Considera-se que este criterio e valido para faixas de triaxialidades medias a 
moderadas. Para triaxialidades baixas ou negativas nao e recomendado aplicar este 
criterio, pois a curva inteira de deformagao de falha real nao pode ser ajustada por 
apenas uma unica fungao exponencial, de forma que na pratica o ajuste de e c para 

o criterio de falha de Johnson-Cook e feito para triaxialidades medias a moderadas. 
Deste modo, este criterio e bastante utilizado na avaliagao de impacto em estruturas 
(Grytten et al., 2009a; Dey et al., 2007), e em estudos relacionados a limite ballstico 
(Grytten et al., 2009b). 


3.3.4 Criterio de falha de Wilkins 


Este criterio e baseado na seguinte variavel cumulativa de dano: 



0 


1 

(l ~a<r m Y 


(2 -Afde p 


(25) 
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onde A = max 


V 5 1 s u 


.s’,, 5 2 e s 3 sao as tensoes desviadoras principals, o m e a 


tensao hidrostatica e os parametros a, A e n sao parametros do criterio de falha. 
Considera-se que um elemento falhou quando D > D cnt , e quando estas duas 
condigoes forem satisfeitas simultaneamente: o m < p cr e £ p >0, onde p cr e uma 
tensao hidrostatica critica, dada por: 


Per 



(26) 


3.3.5 Criterio de falha de Cockcroft-Latham 


Neste criterio a falha ocorre quando a seguinte integral atinge um valor crftico 


D =\ 



de 

p 


(27) 


onde (<j,) e uma fungao que e igual a cr,, se cr, >0, e igual a zero se cr, <0. Este 

criterio de falha e impreciso, pois nao representa com eficacia a falha em altas 
triaxialidades (Teng e Wierzbicki, 2006) alem do valor crftico D ml variar bastante 
dependendo do teste utilizado para a calibragao (para calibrar este criterio e 
necessario apenas um teste experimental). Para alumfnio 2024-T351 o valor D crit 

varia de 0,058 ate 0,485 segundo Teng e Wierzbicki (2006). 

Alem disso, o criterio de Cockcroft-Latham (Cockcroft e Latham, 1968) foi 
estudado por Borvik em simulagoes de limite balfstico (Borvik et al., 2009; 2010) e 
perfuragao de placas (Borvik et al., 2011), resultando em boa correlagao entre 
resultados experimentais e numericos. 
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3.3.6 Criterio de falha de Bao-Wierzbicki 


Este criterio de falha possibilita uma melhor representagao da deformagao de 
falha em fungao do fator de triaxialidade, pois a fungao e dividida em subequagoes 
de acordo com a faixa de triaxialidade. Este criterio e baseado na seguinte variavel 
cumulativa de dano: 


D =\ 



(28) 


onde a fungao de deformagao de falha e c para o alumfnio 2024-T351 e dada por: 


£, =< 
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0 < % < 0,4 
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a 
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(29) 


de forma que o m e a tensao hidrostatica e a e a tensao equivalente. Uma 
particularidade deste criterio de falha e que para triaxialidades abaixo de -1/3 nao 
ha contribuigao para a variavel cumulativa de dano. Assim, fica implementado neste 
criterio a observagao experimental de que para triaxialidades abaixo de -1/3 nao e 
possfvel a ocorrencia de falha (Bao e Wierzbicki, 2005). O significado ffsico dessa 
afirmagao e de que uma trinca nunca podera surgir em urn material submetido a 
compressao hidrostatica com a Jo < -1/3. 
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3.3.7 Criterio de falha de Xue-Wierzbicki 


Este criterio de falha leva em consideragao na sua formulagao o angulo de 
Lode, alem do fator de triaxialidade, o que aumenta consideravelmente sua 
acuracidade. A fungao de dano cumulativo para este criterio tem o formato: 



(30) 


tal que e c = C t e Cin ~{c t e Cl!1 -C 3 e C4 ' ? )(l-^ 1/ ")" onde tj e o fator de triaxialidade e £ e 
o parametro de Lode, dado por: 


£ 


27 ■S’jS'j.S'j 

2 (T 3 


(31) 


Os valores C n sao parametros do criterio de falha, e n e o coeficiente de 
encruamento do material. 


3.4 Modelos micro-mecanicos 


Nos modelos micromecanicos, o material e considerado como sendo urn meio 
poroso, visto que observagoes em escala micrometrica indicam que o mecanismo de 
acumulo de dano nos materiais ducteis envolve a nucleagao, crescimento e 
coalescencia de vazios. Varios trabalhos cientfficos abordaram o tema sobre 
nucleagao-crescimento-coalescencia de vazios sujeitos a varios estados de tensao. 
Nesses modelos, o vazio e modelado como uma esfera ou elipsoide rodeado por 
material sem danejamento, o qual obedece as leis da mecanica do contfnuo. Deste 
modo, a resposta global da porgao de material contendo o vazio e considerada como 
sendo a resposta mecanica macroscopica do material. Devido a enorme combinagao 
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de formatos e orientagoes possiveis para os vazios, sao necessarias simplificagoes 
para tornar o problema matematicamente soluvel. Os trabalhos pioneiros cabem a 
McClintock (1968) e Rice e Tracey (1969) que analisaram o comportamento de um 
vazio cilmdrico ou esferico submetido a tensao uniforme. Nesses modelos, o 
material e considerado como um meio poroso, e a falha ocorre quando e atingido um 
valor crftico para a fragao volumetrica de vazios. Gurson (1977) introduziu o 
acoplamento entre a influencia dos microvazios e a plasticidade do material, e 
Tvergaard (1981) e Tvergaard e Needleman (1984) estenderam o modelo de Gurson 
a fim de incluir a aceleragao do crescimento de vazios. 


3.4.1 Modelo de falha de McClintock 


Foi McClintock (1968) que realizou uma primeira analise teorica do 
crescimento de vazios. O material foi considerado como contendo vazios cilfndricos 
de segao transversal eifptica, de acordo com a Figura 16. 


3 


Figura 16 . Modelo de crescimento de vazios de McClintock. 
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McClintock analisou uma longa cavidade cilfndrica em um material tracionado 
no eixo da cavidade, enquanto o mesmo era submetido a tensoes transversais. O 
crescimento do vazio cilmdrico era governado por: 


d log 


Ssenh 


Sc 


r,inf 


y ^r ,inf &Z .inf ) 


inf + ^r, inf 


(32) 
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onde ze a diregao axial, re a diregao radial, /'nfdenota valores distantes, £ 0 e o raio 

original do cilindro e b e o raio instantaneo do mesmo. O modelo de dano de 
McClintock pode ser expresso por: 


dD 


log F : 


zb 


s 


2(1 -n) 


senh\ 


S{\-n)(o a -(J h ) 


+ - 


3 


de. 


(33) 


onde log^{ denota a maxima deformagao da celula contendo o microvazio, cr a e a b 
sao as tensoes transversais, o e a tensao equivalente, e p e a deformagao plastica 

equivalente, n e o coeficiente de encruamento e D e a variavel cumulativa de dano. 
A falha ocorre quando D atinge o valor unitario. 


3.4.2 Modelo de Rice e Tracey 


Rice e Tracey (1969) modelaram o crescimento de urn vazio esferico em urn 
material puramente elastico. Urn vazio esferico sujeito a tensao uniaxial nao 
somente cresce em tamanho, mas tambem muda de forma. Rice e Tracey 
mostraram que quanto maior o fator de triaxialidade, mais significativa e a mudanga 
de tamanho do vazio em comparagao com a sua mudanga de forma. Essa 
conclusao esta de acordo com a afirmagao de que a tensao hidrostatica tern por 
efeito acelerar o processo de nucleagao-crescimento-coalescencia de vazios. Rice e 
Tracey obtiveram a taxa de crescimento do raio dos vazios como sendo: 


— = 0,283 exp| 

Rn 


r 3(7,,, ^ 


2(7 


esc J 


(34) 


onde R 0 e o raio original do vazio, Re o seu valor atual, <7 m e a tensao hidrostatica 
medida a distancia e <r esc e a tensao de escoamento. Segundo os autores, para urn 
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material sofrendo encruamento sob tensao proporcional, a deformagao de falha 
pode ser expressa por: 


c/cxp 


— b- 


eq y 


(35) 


onde a e b sao dois parametros do modelo. Neste modelo, Rice e Tracey nao 
consideraram a interagao entre vazios adjacentes. 


3.4.3 Modelo de Gurson-Tvergaard-Needleman 


O modelo de plasticidade com efeito de porosidade e chamado de modelo de 
Gurson-Tvergaard-Needleman. Ele e baseado no modelo constitutivo poroso de 
plasticidade, que e dependente da tensao hidrostatica e foi desenvolvido por Gurson 
(1977) considerando nucleagao e crescimento de vazios. Considera-se aqui que os 
vazios estao dispersos em urn meio continuo incompressivel. O parametro 
denominado fragao volumetrica de vazios ^especifica o imcio do surgimento de uma 
trinca, quando o parametro atinge urn valor critico. 

A fungao de escoamento proposta por Gurson e dada por: 


<t> = + 2/ cosh 


f ^ „ \ 

3 <y,„ 


2 (7. 


- 1 -/, =0 


(36) 


onde <y e a tensao equivalente de von Mises, a m e a tensao hidrostatica e a esc e a 

tensao de escoamento do material. Posteriormente foi realizada uma modificagao no 
modelo de Gurson de forma a introduzir dois novos parametros qi e Assim, a 
fungao de escoamento torna-se: 


<T = 


rf 2 


+ 2 qj v cosh 


V 2 °escj 


V 


1 -(«,/„ ) 2 =0 


(37) 
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de acordo com a equagao (37), o material perde capacidade de transmitir tensao 

quando os vazios apresentam grande crescimento. A fungao f (/v) foi criada por 
Tvergaard e Needleman (1984) para substituir f v na equagao anterior, a fim de 
modelar a perda gradativa de resistencia do material: 

antes de qualquer carregamento 
para f v < f vc 

Para f v > f vc 

(38) 

onde f v0 e a fragao volumetrica de vazios inicial, f vc e o valor crftico da fragao 
volumetrica de vazios, a partir do qual se forma uma trinca (e a capacidade de 
transmitir carga cai rapidamente) e f vF e a fragao volumetrica de falha na qual 

ocorre perda completa da resistencia do material (ruptura completa gerada por 
propagagao de trinca). 

O incremento da variavel de fragao volumetrica de vazios resulta da 
nucleagao e coaiescencia de vazios: 


/'(/,)= 


/v0> 

fv, 

-fvc 

/«+4-^(/v-/J 


1 


fvF fvc 


d.f\, = dj\ 


v _ nucleagao 


+ df\ 


v _ coaiescencia 


(39) 


Gurson propos que a taxa do crescimento da fragao volumetrica de vazios 
devido a coaiescencia e dada por: 


df\, 


coaiescencia 




(40) 


p 

onde e a taxa de deformagao plastica volumetrica. Chu e Needleman (1980) 
consideraram uma abordagem estatfstica para a nucleagao de vazios. Os autores 
propuseram uma distribuigao normal das nucleagoes de vazios com relagao a 
deformagao plastica: 
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f — A £ 

J v _ nucleagao N*’' p 


onde s p e a deformagao plastica na matriz cristalina e 


F 

A - N 


S n ^27T 


exp- 


f \2 

£ P ~£,, 


\ J 


(41) 


(42) 


onde Sn e s n sao o desvio padrao e a media da distribuigao de deformagoes 
plasticas e F N e a fragao volumetrica total que pode ser nucleada. 

Apesar do sucesso dos modelos de nucleagao-crescimento-coalescencia de 
vazios na predigao de falha em varias aplicagoes, ainda existem empecilhos com 
estes tipos de modelo: 

- nao e tratado o mecanismo de cisalhamento de vazios. O modelo de Gurson falha 
na simulagao de urn simples teste de cisalhamento. O modelo original de Gurson foi 
concebido para representar bem a faixa de alta triaxialidade, onde a pressao 
influencia diretamente no fenomeno de falha. Assim, recomenda-se cuidado na 
utilizagao do modelo de Gurson em baixa triaxialidade. 

- no modelo de Gurson o danejamento e restrito a fragao volumetrica de vazios. No 
entanto, o processo de danejamento no material e mais complexo do que a evolugao 
da distribuigao de vazios. Por exemplo, uma trinca sem volume na estrutura 
cristalina do material pode introduzir danejamento ao material ao fazer o mesmo 
perder resistencia mecanica. 

- trabalhos posteriores adicionaram novas formulas ao modelo de Gurson, no 
entanto, o modelo ainda nao se mostra pratico em termos de facilidade de 
calibragao. De fato, com os desenvolvimentos posteriores do modelo, buscando 
aumentar a sua aplicabilidade, chegou-se a urn total de dez parametros a serem 
determinados, sendo que os mesmos estao fortemente acoplados. Por essa razao 
encontram-se poucos casos de aplicagao desse modelo na industria. 
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3.5 Modelo de dano de Lemaitre (mecanica do dano contfnuo) 

Modelos de dano que detectam a deterioragao do material desde o infcio do 
escoamento inserem uma variavel de dano D a lei constitutive, que varia de 0 a 1, 
onde D = 0 representa material fntegro e D = 1 representa perda total da capacidade 
de carga. O dano resulta na perda de rigidez, possibilitando que a falha se inicie na 
regiao danificada. 

A diferenga dos modelos de mecanica do dano contfnuo em comparagao com 
os modelos micromecanicos esta no fato de o tratamento do crescimento e interagao 
de vazios do primeiro ser tratado fenomenologicamente. Deste modo, parametros 
macroscopicos sao utilizados na descrigao da resposta mecanica do material. 

No modelo de dano de Lemaitre, a lei de elasticidade com dano acoplado e 
dada por: 


a = {l-D )D e :f 


(43) 


onde De a variavel de dano, D 6 e o tensor de elasticidade isotropico e s e e o tensor 
de deformagao elastico. A lei de evolugao de dano e dada por: 



(44) 


onde S e s sao parametros do material, e Y e o multiplicador plastico, urn fator 


relacionado ao incremento de plasticidade, calculado a cada ciclo de integragao no 
tempo (ver capftulo 5). A variavel -Y, que representa a forga termodinamica 
associada com o dano, e igual a: 



( 45 ) 
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onde crea tensao equivalente de von Mises, pea pressao hidrostatica, Geo 

modulo cisalhante e Keo modulo volumetrico. 

Como exemplo adicional, a variavel de dano da mecanica do dano continuo 
tambem pode ser introduzida em modelos constitutivos, de forma acoplada, como 
por exemplo no modelo de plasticidade de Johnson e Cook (1985): 


5 = (l-P b D\A+Be- 


1 + C In — 

£ n 


1 - 


T — T 
1 1 o 


T f 

V / 


l o J 


(46) 


onde D e a variavel de dano cumulativo do criterio de falha de Johnson-Cook e Pb e 
urn parametro de material. 
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4 ANALISE EXPERIMENTAL 


Neste trabalho propoe-se o projeto da geometria de um corpo de prova 
especial que falhe em diferentes rnveis de triaxialidade, especificamente alta 
triaxialidade (falha por nucleagao, crescimento e propagagao de vazios) e 
triaxialidade proxima de zero (falha nas bandas de cisalhamento). 


4.1 Geometria inicial 


A geometria inicial do corpo de prova foi proposta por Alves e Driemeier 
(2010), conforme Figuras 17 e 18. 



Figura 17 . Primeira geometria do especime Bifailure, e corpo de prova na maquina de tragao. 
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Figura 18. Medidas da primeira geometria do especime Bifailure 


Para a primeira geometria do Bifailure foi testado o material alummio 
aeronautico, dispomvel em chapas de 6,35mm de espessura, com caracteristicas 
dadas naTabela 1. 


Tabela 1 - Propriedades do material alummio aeronautico 


Densidade 

Modulo de 

Coef. de Poisson 

E 

E 

o 

H— 1 

elasticidade (MPa) 


2,78x10 -9 

67100 

0,33 


A Figura 19 ilustra a sequencia de instantes do ensaio da primeira geometria 
do especime. Percebe-se que, apos a primeira ruptura no entalhe menor (regiao de 
alta triaxialidade), houve grande rotagao da regiao central do especime, afastando a 
falha na regiao central da condigao desejavel de cisalhamento puro. 
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Figura 19. Resultados experimentais para a primeira geometria do Bifailure, em aluminio 


A Figura 20 mostra claramente a plastificagao consideravel de regioes fora 
da regiao de ruptura, ocasionadas, principalmente, pelo desequilfbrio do especime e 
rotagao do mesmo apos a primeira falha. 



Figura 20. A grande deformagao das regioes destacadas dificulta a obtengao de cisalhamento puro 

na regiao central do especime. 
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4.2 Geometria proposta 


Como solugao para o problema de rotagao do especime, foi proposta uma 
nova geometria para o especime Bifailure, com dois entalhes menores, conforme 
iiustra a Figura 21. A geometria do especime foi projetada para que a deformagao de 
falha dos entalhes menores fosse pequena e, portanto, nao houvesse rotagao 
significativa do especime. As medidas do especime sao apresentadas nas Figuras 
22 e 23. Este especime com as dimensoes mostradas a seguir foi o especime 
utilizado neste trabalho. 



Figura 21. Segunda proposta do especime Bifailure 
145 
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Figura 23. Em a), distanciamento do centro dos furos dos entalhes menores em relagao a parede 
lateral do especime; em b), indicagao de raio e espagamento maximo da curvatura dos entalhes em 

relagao a face superior. 

Os ensaios de tragao foram realizados na maquina de tragao Instron do 
Laboratorio de Impacto do GMSIE, conforme ilustrado na Figura 24. 



Figura 24. Corpo de prova montado na maquina de tragao 


Para os testes foi utiiizado um extensometro para medigao de deformagao, 
ilustrado na Figura 25. 
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Figura 25. Corpo de prova Bifailure montado na maquina de tragao com extensometro 


4.2.1 Caracterizagao do material Alummio 2024-T351 


O material utilizado neste trabalho foi a liga de alummio 2024-T351, de forma 
que os corpos de prova foram fabricados a partir de uma chapa com 9,5 mm de 
espessura do referido material. As propriedades desse material sao listadas na 
Tabela 2. 


Tabela 2 - Propriedades do material Alummio 2024-T351 


Densidade 

Modulo de 

Coef. de Poisson 

(ton/mm 3 ) 

elasticidade (MPa) 


2,7x10 -9 

74660 

0,3 


Foram realizados testes de tragao uniaxial de corpo de prova padrao para 
esse tipo de ensaio, com dimensoes dadas pela Figura 26. 
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32 <5 R = 12,5 
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Figura 26. Medidas do corpo de prova para ensaio de tragao 

O ensaio de tragao padrao consiste em tracionar o corpo de prova ate a 
ruptura. O corpo de prova utiiizado e ilustrado na Figura 27. 



Figura 27. Corpo de prova padrao utiiizado para ensaio de tragao 


A curva tensao-deformagao de engenharia obtida para o alurmnio 2024-T351 
e ilustrada na Figura 28. Verifica-se na curva o instante de ruptura do especime, 
evidenciado pela queda da magnitude da curva no final do teste. 

A opgao pelo estudo da liga de alurmnio 2024-T351 foi pelo fato de os 
parametros que definem as propriedades do material estarem dispomveis na 
literatura cientffica. Trabalhos como o de Wierzbicki et al. (2005) e Malcher et al. 
(2012) apresentam urn extensivo estudo sobre modelos de falha apiicados ao 
alurmnio 2024-T351 e portanto fornecem tabelas com os parametros dos modelos 
de falha utilizados neste trabalho. 
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Figura 28. Curva tensao-deformagao de engenharia para aluminio 2024-T351 



Figura 29. Resultado experimental dos testes de tragao padrao 


Posteriormente foram realizados testes de tragao do especime Bifailure 
(segunda proposta), obtendo-se o resultado mostrado na Figura 30. Na mesma 
figura aparece tambem o indicativo da regiao de falha correspondente a cada queda 
da curva de forga do grafico. Desse modo, obtem-se que a primeira queda da curva 
de forga esta ligada a ruptura dos entalhes laterais menores, enquanto que a queda 
final da curva de forga esta ligada a ruptura do entalhe central. 
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Obteve-se ruptura do especime Bifailure sem grande rotagao da regiao 
central, de acordo com a Figura 31, o que favoreceu a obtengao de medidas mais 
precisas com o extensometro. 



Figura 31. Regiao central do especime Bifailure apos o ensaio de tragao 


Em smtese, as propriedades do alummio 2024-T351 utilizadas neste trabalho 


sao listadas na Tabela 3. 
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Tabela 3 - Propriedades para aluminio 2024-T351 


Descripao 

Simbolo 

Valor 

Referenda 

Densidade 

P 

2700 kg/m 3 

Teng et al. 2005 

Modulo elastico 

E 

74660 MPa 

Teng et al. 2005 

Coef. Poisson 

V 

0,3 

Teng et al. 2005 

Parametro de JC 

A 

352 MPa 

Teng et al. 2005 

Parametro de JC 

B 

440 MPa 

Teng et al. 2005 

Parametro de JC 

n 

0,42 

Teng et al. 2005 

Max. def. plast. equiv. de falha 


0,21 

Teng e Wierzbicki, 2006 

Max. tensao cisalhante de falha 

T 

max 

380 MPa 

Teng e Wierzbicki, 2006 

Parametro de falha de JC 

A 

0,13 

Teng et al. 2005 

Parametro de falha de JC 

d 2 

0,13 

Teng et al. 2005 

Parametro de falha de JC 

D, 

-1,5 

Teng et al. 2005 

Parametro de falha de Wilkins 

a 

0,0012 MPa 1 

Wierzbicki et al. 2005 

Parametro de falha de Wilkins 

X 

2,15 

Wierzbicki et al. 2005 

Parametro de falha de Wilkins 

P 

2,18 

Wierzbicki et al. 2005 

Parametro de falha de Wilkins 

D c 

0,93 

Wierzbicki et al. 2005 

Valor crftico superior de falha de 

Dent 

0,485 

Teng e Wierzbicki, 2006 

Cockcroft-Latham 




Valor crftico inferior de falha de 

Dcrit 

0,058 

Teng e Wierzbicki, 2006 

Cockcroft-Latham 




Parametro de falha de Xue-Wierzbicki 

Ct 

0,87 

Wierzbicki et al. 2005 

Parametro de falha de Xue-Wierzbicki 

c 2 

1,77 

Wierzbicki et al. 2005 

Parametro de falha de Xue-Wierzbicki 

c 3 

0,21 

Wierzbicki et al. 2005 

Parametro de falha de Xue-Wierzbicki 

C4 

0,01 

Wierzbicki et al. 2005 

Coeficiente de encruamento 

n 

0,153 

Wierzbicki et al. 2005 
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5 IMPLEMENTAQAO NUMERICA 


5.1 Implementagao de subrotina 


A fim de se obter maior controle sobre o algoritmo numerico de solugao e para 
ser possivel a utilizagao de criterios de falha que nao estao originalmente 
dispomveis no programa comercial de elementos finitos, foi desenvoivida uma 
subrotina de material para o software LS-Dyna®. A subrotina consiste na 
implementagao do modelo constitutivo de material e dos criterios de falha que foram 
tratados neste trabalho. O modelo de material que foi implementado foi o de von 
Mises, para lidar com elementos solidos. 

A subrotina e escrita na linguagem de programagao Fortran, a qual e tambem 
utilizada originalmente na compilagao do programa de elementos finitos LS-Dyna®. 

A subrotina completa esta disponfvel no Apendice B, sendo que aqui serao 
explicados os pontos mais importantes sobre o seu funcionamento. 

O modelo de plasticidade de von Mises e composto por: 

1) Uma lei linear elastica: 


rx = D e : /; 


(47) 


onde Deo tensor de elasticidade isotropico e s e e o tensor de deformagoes. 


2) Uma fungao de escoamento na forma: 

c y ) = fiJJs) -o,= a(s)- cr y 


(48) 


onde s e o tensor desviador, o e a tensao equivalente e <x v e a tensao de 
escoamento, que e tomada como uma fungao da deformagao plastica acumulada, 
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ou fungao do modelo constitutivo, de forma que neste trabalho e igual a 
< 7 y = A + Be" (modelo constitutivo de von Mises). 

3) Uma lei de escoamento: 



(49) 


tal que y e chamado de multiplicador plastico. 

4) E uma lei de encruamento, com a equagao de evolugao para a variavel interna de 
encruamento: 



(50) 


O processo de retorno aplicado ao estado de previsao elastica descrito na 
proxima segao e necessario para se determinar quanto do incremento de 
deformagao total que ocorre a cada passo de integragao e de natureza elastica e 
quanto e de natureza plastica. O fluxograma da Figura 32 ilustra o que se explicou 
aqui. Nesse fluxograma, o fndice n indica o passo atual e o fndice n+1 indica o 
proximo passo no processo de integragao. 
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Figura 32. Fluxograma representando o funcionamento da solugao numerica de elementos finitos 


5.2 O procedimento de previsao elastica / algoritmo de retorno 


Dado o incremento de deformagao total: 

A e n+l = e n+1 -e„ (51) 


correspondente a um pseudo incremento de tempo [t n ,t n+1 ] e dadas as variaveis de 
estado {e®,} em t n , a fim de ser possfvel determinar quanto do incremento de 
deformagao total e elastico e quanto e plastico considera-se inicialmente que esse 
incremento de deformagao total e elastico. Desse modo, a deformagao elastica de 
tentativa e a deformagao plastica acumulada de tentativa sao dadas por: 


c e,tent 

b n+l 


= + As 


(52) 


p p,tent _ pP 
C n+1 C n 


( 53 ) 
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A correspondente tensao de tentativa e caiculada por: 


<n7, = d 


\s 


n+l 


(54) 


Ou, de forma equivalente, decompondo nas partes hidrostatica e desviadora: 


tent 
s n+l 


= 2Gs 


e .tent 
d.n+1 


P 


tent 

n+l 


= k>: 


e.tent 

v.n+1 


(55) 

(56) 


onde sep denotam, respectivamente, os tensores desviador e hidrostatico, G e K 
sao respectivamente os modulos cisalhante e volumetrico e os indices dev 
representam, respectivamente, os componentes desviador e volumetrico. 

Depois de calcular o estado elastico de tentativa, o proximo passo no 

— tent 

algoritmo e verificar se °n+i se situa dentro ou fora da superffcie de escoamento de 
tentativa: 

— tent 

• Se °n+ise situar dentro da superffcie de escoamento de tentativa, isto e, se: 


0(<T f ^,a y!M ) =o-fc)-^ >n (57) 

entao o processo no intervalo [t n ,t n+1 ] e puramente elastico e o proprio estado 
elastico e a solugao do problema de integragao. Desse modo, se atualiza: 


a e — p e J ent 

t'n+l ~ 1 


_, _ — tent 

^ n+l ~ ^ n+l 


£ p =£ p 
c n +1 c n 


<J ^=<7 

y,n +1 y,n 


(58) 


• Caso contrario, o processo e elasto-plastico no intervalo [t n ,t n+1 ] e o algoritmo de 
retorno descrito a seguir precisa ser aplicado. A Figura 33 ilustra este caso, 
mostrando o processo de previsao elastica e algoritmo de retorno a superffcie de 
escoamento. 
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Figura 33. Representagao esquematica da predigao elastica / algoritmo de retorno para o modelo de 

von Mises. Extraido de Souza Neto et al. (2008). 


Deve-se destacar que o algoritmo de retorno consiste em resolver o sistema 
de equagoes nao lineares a seguir: 


Cl = e;,T - Ay, 


3 


2 \\s. 


:+l 


E n+\ =e;+Ay 

V 3 A(0-°,fe«)= o 


Que precisa ser resolvido para e‘ +1 . C+i e Ay, onde: 


S n+1 S fj+l( e n+l) “^^d.n+1 


(59) 


(60) 


Depois de resolver o sistema de equagoes (59), o tensor de deformagoes plasticas 
pode ser atualizado de acordo com a seguinte formula: 


< + 1 =e'+AyJ- 


3 s„ 


2 s 


M+l 


(61) 
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O sistema de equagoes (59) pode ser simplificado em uma unica equagao 
nao-linear em fungao do multiplicador plastico Ay. Esta redugao no numero de 
equagoes e fundamental para tornar o processo de atualizagao de tensoes mais 
eficiente computacionalmente. Pela simplificagao do sistema de equagoes (59), 
nota-se que a plasticidade afeta apenas a parte desviadora e nao a volumetrica dos 
tensores. A separagao em partes desviadora e volumetrica das deformagoes 
fornece: 


c e _ ~e,tent 

^v,n+l — ^v,n+l 


Ay J- 


3 s„ 


2 s„ 


De forma equivalente, em termos de tensoes temos: 


Pn"l 


•V„_, =*n -Ay 2 GJ^^L 


3 s„ 


2 s„ 


(62) 

(63) 


(64) 

(65) 


Assim, o processo de retorno afeta somente o componente desviador de 
tensao. A componente hidrostatica, p n+l , tern seu valor calculado na fase de 

previsao elastica e pode, portanto, ser eliminada do sistema de equagoes. Uma 
maior simplificagao e alcangada rearranjando a equagao (65): 


f 

1 + 

V 


3 Ay2 G 


y n +1 


n+ 11 | 


tent 

S n+1 


( 66 ) 


ou seja, as tensoes desviadora de tentativa e atualizada sao colineares. Isso implica 
que: 


S 


n +1 


tent 

*n +1 



( 67 ) 
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A substituigao da equagao (67) na equagao (65) conduz a uma formula mais simples 
para a atualizagao de tensoes desviadoras: 



«+1 


1- 


Ay3 G 

_ tent 

o n +i y 


'w+1 


( 68 ) 


— tent = /sjr Arenr ) 

onde <T,1+1 v 2 V »+i/ e a tensao elastica de tentativa de von Mises. Observa-se 

tent 

que sendo urn tensor constante no processo de retorno, a tensao desviadora 
e uma fungao linear de Ay somente na equagao de atualizagao (68). Temos que 
a equagao (68) indica que a tensao desviadora atualizada e obtida aplicando-se o 
fator de escala l-A/3G/a' n e * a tensao desviadora de tentativa. 

Finalmente, com a substituigao das equagoes (68) e (59) 2 na equagao (59) 3 , 
obtem-se que o sistema de equagoes (59) e reduzido a seguinte equagao escalar (e 
nao linear) que e fungao do multiplicador plastico Ay. 


*(A /) = <?“ -3GA y-ojpl +Ar)=0 (69) 

A equagao (69) e resolvida pelo metodo de Newton-Raphson, e com a solugao Ay, 
as variaveis de estado sao atualizadas de acordo com: 


S n + 1 = 


1- 


Ay3 G 

_ tent 

& n+1 7 


n+1 


_ „ , tent -r 

CT «+1 — S n +1 ■*" Pn+ 17 

[Vve l - ^ f „ . f „e.tent 

- [D J -<* n+1 — ^^ s n +l + ^ £ v.n+l 


8 n + i = e „+Ay 


(70) 


O procedimento de estado elastico de tentativa e algoritmo de retorno esta 
resumido nos esquemas das Figuras 34 e 35 ao final deste capftulo. 

Os criterios de falha implementados na subrotina sao os mesmos descritos na 
segao 3.3, ou criterios de falha empfricos baseados na mecanica do dano. Vale 
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ressaltar que os criterios de falha implementados produzem leves diferengas de 
resultados em comparagao com os mesmos criterios ja dispomveis no programa de 
elementos finitos, como o criterio da maxima deformagao plastica equivalente, ou 
criterio de falha de Johnson-Cook. Acredita-se que essas diferengas ocorrem pelo 
criterio de ativagao de falha nos elementos. Na subrotina, conhece-se que a falha e 
ativada quando o valor crftico de falha e atingido em pelo menos urn ponto de 
integragao do elemento. 
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Algoritmo de previsao elastica e retorno 


(i) Previsao elastica. Dado Ae e as variaveis de estado em t n , avaliar o 
estado elastico de tentativa: 



(ii) Verificar ocorrencia de plasticidade: 

SE 

ENTAO ajuste (.) n+1 := (.)*" e SAIA da subrotina 

(iii) Procedimento de retorno. Resolver a equagao: 

= a';:; -3CAy-<j y & +Ay)=0 

para Ay usando o metodo de Newton-Raphson (ir para esquema da Figura 
35) e atualizar as variaveis de estado: 

n ~ n ,aU ■ s fl _ .AjjgU, 

rn+\ m Fn+ 1 » 3 n+l • 1 —tent 3 h+1 

< °”«+1 / 

: = *„ + l + Pn + J 



(iv) SAIR 


Figura 34. Esquema do algoritmo de previsao elastica e retorno. Extrafdo de Souza Neto et al. 

( 2008 ). 
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Algoritmo de Newton-Raphson 

(i) Ajuste um valor inicial para Ay 

A/ 0) := 0 

e o correspondente residual (valor da fungao de escoamento): 

(ii) Executar iteragao de Newton-Raphson 

ch 7 V 

H := —^ (inclinagao) 

de p A 

e„ p +Af 

d :=-= -3G -// (derivada residual) 

Ay:= Ay - (novo valor para Ay) 

d 

(iii) Verificar convergence 

$ := a 1 ::; - 3GAy-<r y (e p + Ay) 

SE d <s tol ENTAO RETORNAR ao esquema da Figura 34 

(iv) IR PARA (ii) 

Figura 35. Esquema do algoritmo de Newton-Raphson. Extraido de Souza Neto et al. (2008). 
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6 ESTUDO DE CASO DA APLICAQAO DOS CRITERIOS DE FALHA 


Conforme explicado anteriormente, a fim de simular numericamente os 
ensaios experimentais, foi utilizado o programa comercial de elementos finitos LS- 
Dyna i® e o modelo constitutivo de von Mises: 


g = A + Bs n p 


( 71 ) 


onde A, B e n sao constantes do material que precisam ser determinadas 
experimentalmente. Neste trabalho esses parametros foram extrafdos do modelo 
constitutivo de Johnson-Cook, sem considerar os efeitos da taxa de deformagao e 
temperatura, a partir do trabalho de Teng et al. (2005). 


Tabela 4 - Parametros do modelo de von Mises para alummio 2024-T351 


A (MPa) 

B (MPa) 

n 

352 

440 

0,42 


Inicialmente foi realizada a simulagao do teste de tragao com especime 
comum, utilizando a malha de elementos finitos mostrada na Figura 36. 



Figura 36. Malha de elementos finitos utilizada na simulagao do teste de tragao comum 


A simulagao do teste de tragao comum foi feita a fim de verificar a correlagao 
entre as curvas tensao-deformagao experimental e numerica. Nesta comparagao foi 
inclufda tambem a curva do alummio 2024-T351 utilizada por Wierzbicki et al. 
(2005). 
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Deformagao verdadeira 

Figura 37. Curva tensao-deformagao para teste de tragao comum, com a curva extrafda da simulagao 
e curva do material utilizado por Wierzbicki et al. (2005) 

Verifica-se na Figura 37 boa correlagao entre resultados experimentais e 
numericos, com relagao ao comportamento plastico do material. 

Como ocorreu a validagao dos parametros de material adotados, os mesmos 
foram utilizados na simulagao do especime Bifailure, no teste de diferentes criterios 
de falha. A malha do especime Bifailure e apresentada nas Figuras 38 e 39. Os 
pinos que tracionam o especime foram modelados como superficies analfticas 
(cilindros), sendo aplicado urn deslocamento imposto no pino de uma extremidade 
enquanto que o outro pino e mantido fixo. Os valores de deformagao de engenharia 
obtidos sao tornados a partir da distancia entre dois nos que corresponderiam aos 
pontos de apoio das pingas do extensometro com o corpo de prova. Por 
simplificagao, os entalhes que rompem a alta triaxialidade serao chamados de 1 e 2, 
e o entalhe central, que rompe a triaxialidade proxima a zero sera chamado de 
entalhe 3. 



Figura 38. Malha do corpo de prova Bifailure , composta de 152220 elementos 
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Figura 39. Detalhe da malha do especime Bifailure, na regiao dos entalhes 


Atraves de simulagao foram medidas as triaxialidades para os entalhes 1,2 e 
3, que estao ilustradas na Figura 40. O valor medio de triaxialidade obtido para os 
entalhes 1 e 2 e de 0,629 e para o entalhe 3 e de 0,016. 



Figura 40. Variagao da triaxialidade para os entalhes do especime Bifailure 


A seguir sao discutidos os resultados numericos dos criterios de falha 
implementados na subrotina de elementos finitos. 
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6.1 Criterio de falha da maxima deformagao plastica equivalente 


Para este criterio foi utilizado o valor de deformagao de falha de 0,21, extrafdo 
do trabalho de Teng e Wierzbicki (2006). O resultado da simulagao e apresentado na 
Figura 41. 



Deformagao de engenharia (%) 

Figura 41. Curvas experimental e numerica da simulagao com criterio de falha da maxima 

deformagao plastica equivalente 

Verifica-se correlagao razoavel entre resultados experimentais e numericos, 
apesar do criterio embasado apenas em deformagao ser relativamente simples. 
Apos a ruptura dos entalhes 1 e 2, verifica-se na curva numerica urn comportamento 
oscilatorio do especime. Isso ocorre pelo fato de nao haver amortecimento nas 
simulagoes. 



Figura 42. Resultado da simulagao com criterio de falha da maxima deformagao plastica equivalente 
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Verifica-se que a acuracidade deste criterio depende do valor utilizado para a 
calibragao do mesmo. Pelo fato de ser necessario apenas um teste para a 
calibragao do criterio da maxima deformagao plastica equivalente, em Wierzbicki et 
al. (2005) foi adotado um valor que produzisse uma melhor correlagao com os testes 
experimentais analisados nesse trabalho. O valor utilizado para a calibragao nao foi 
baseado no teste de tragao uniaxial, o qual produziria um valor crftico muito maior 
(£ f =0,45). 


6.2 Criterio de falha da maxima tensao cisalhante 


Para este criterio foi utilizado o valor de tensao cisalhante de falha de 380 
MPa, extrafdo do trabalho de Teng e Wierzbicki (2006). O resultado da simulagao e 
apresentado na Figura 43. Verifica-se que o criterio da maxima tensao cisalhante 
superestimou a ruptura do especime, resultando em grande rotagao da parte central 
do mesmo, evidenciada na Figura 44, que ilustra o momento imediatamente apos a 
falha. 



Figura 43. Curvas experimental e numerica da simulagao com criterio de falha da maxima tensao 

cisalhante 
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Figura 44. Resultado da simulagao com criterio de falha da maxima tensao cisalhante 


Este criterio precisa apenas de um teste experimental para ser calibrado, por 
essa razao o valor utilizado para a calibragao deve ser o mais representative 
possivel da ruptura dos especimes ensaiados. Uma implementagao do criterio da 
maxima tensao cisalhante foi feita no espago das deformagoes principals em 
Wierzbicki et al. (2005), chegando-se a bons resultados. 


6.3 Criterio de falha de Johnson-Cook (JC) 


Para este criterio, a deformagao de falha tern o seguinte formato: 

£ c = d i + D i exp {D 3 tj) ( 72 ) 

Desta forma, neste trabalho foram adotados os valores de Teng et al. (2005), 
que aparecem na Tabela 5. A relagao fator de triaxialidade x deformagao de falha 
esta mostrada na Figura 45. 

Tabela 5 - Parametros do criterio de falha de Johnson-Cook para alurmnio 2024-T351 


Di 

d 2 

d 3 

0,13 

0,13 

-1,5 
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Fator de triaxialidade V 

Figura 45. Lugar geometrico de falha para o criterio de falha de Johnson-Cook 



Deformagao de engenharia (%) 

Figura 46. Resultado da simulagao com criterio de falha de Johnson-Cook 

De acordo com a Figura 46, verifica-se boa correlagao entre resultado 
experimental e numerico, apesar de o criterio de falha de Johnson-Cook considerar 
apenas o fator de triaxialidade em sua formulagao. Tres testes experimentais sao 
necessarios para a calibragao do criterio de falha de JC. Desse modo, da mesma 
forma que ocorre com os dois criterios anteriores, o criterio de JC tambem necessita 
de uma calibragao que satisfaga a maior parte dos ensaios experimentais 
analisados. Por esse motivo encontra-se dispomvel na literatura mais de urn 
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conjunto de parametros de JC para alummio 2024-T351. A Figura 47 ilustra o 
especime logo apos a falha. 



Figura 47. Resultado da simulagao com criterio de falha de Johnson-Cook 


Um ponto interessante a ser considerado e que nesta simulagao a falha do 
entalhe 3 nao ocorreu exatamente no centra, mas sim de forma assimetrica. Outro 
ponto importante foi uma representagao razoavel da ruptura no entalhe 3. O criterio 
de Johnson-Cook foi originalmente proposto para a representagao de falha em altas 
triaxialidades. Assim, pode-se dizer que a representagao de falha no entalhe 3 foi 
razoavel apesar do fato de a ruptura proxima a cisalhamento puro implicar em altos 
valores de deformagao de falha pelo criterio de JC. 


6.4 Criterio de falha de Wilkins 


O criterio de falha de Wilkins implementado aqui foi calibrado de acordo com 
Wierzbicki et al. (2005). Os parametros utilizados aparecem na Tabela 6. 


Tabela 6 - Parametros do criterio de falha de Wilkins para alurmnio 2024-T351 


a 

A 

A 

D C rit 

1,2 10“ 3 MPa~ l 

2,15 

2,18 

0,93 
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Pelos resultados apresentados, na Figura 48 percebe-se que o criterio de 
falha de Wilkins subestimou a ocorrencia de ruptura em alta triaxialidade (entalhes 
laterais menores) e gerou uma pessima previsao da falha em cisalhamento (entalhe 
central). Apesar de serem necessarios no rmnirno 4 testes experimentais para a sua 
calibragao, este criterio de falha nao apresentou bons resultados no presente 
trabalho. De acordo com Wierzbicki et al. (2005) este criterio produziria bons 
resultados na faixa de triaxialidade -1/3 <tj <1/3 enquanto que subestimaria 
resultados na faixa de triaxialidade 1/3 <r/< 2/3. 



Figura 49. Resultado da simulagao com criterio de falha de Wilkins. 
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6.5 Criterio de falha de Cockcroft Latham 


Este criterio de falha e baseado em uma energia de deformagao de ruptura 
na definigao do valor critico de falha. Desse modo, apenas um teste experimental e 
necessario para a sua calibragao. No entanto, da mesma forma que para os criterios 
baseados na maxima deformagao plastica equivalente e no criterio da maxima 
tensao cisalhante, e necessario calibrar este criterio com o valor que melhor 
representa na media os fenomenos de falha analisados. De acordo com Teng e 
Wierzbicki (2006) os valores crfticos deste criterio para o alurmnio 2024-T351 variam 
de 0,058 ate 0,485. Dessa forma, aqui foram feitas simulagoes com esses dois 
valores extremos. 



Deformagao de engenharia (%) 


Figura 50. Resultado da simulagao com criterio de falha de Cockcroft-Latham e valor critico de 0,058 
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Deformagao de engenharia (%) 

Figura 51. Resultado da simulagao com criterio de falha de Cockcroft-Latham e valor critico de 0,485 

Verifica-se que os valores extremos de calibragao deste criterio nao geraram 
bons resultados para a simulagao do especime Bifailure. A fim de se representar 
melhor a simulagao do Bifailure, um melhor valor critico poderia ser escolhido para a 
simulagao. De acordo com Teng e Wierzbicki (2006), a calibragao deste criterio a 
partir de ensaios de compressao de cilindros oferecem valores consistentes, 
enquanto que a calibragao feita a partir de ensaios de tragao uniaxial de especimes 
com e sem entalhe geram valores que se espalham em uma longa faixa. 



Figura 52. Resultado da simulagao com criterio de falha de Cockcroft-Latham e valor critico de 0,058 

E interessante notar na Figura 52 que a utilizagao de um valor critico baixo 
resultou em duas pequenas trincas, uma de cada lado ao redor do entalhe 3, antes 
da ocorrencia da ruptura. Nas Figuras 52 e 53, percebe-se que houve ruptura nao 
simetrica em relagao ao entalhe 3, sendo que na Figura 54 verifica-se a existencia 
de duas microtrincas de modo que uma delas ira originar a ruptura final do entalhe 3 
(microtrinca do lado vermelho). 
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Figura 53. Resultado da simulagao com criterio de falha de Cockcroft-Latham e valor crftico de 0,485 



Figura 54. Trincas ao redor do entalhe 3, para criterio de falha de Cockcroft-Latham com valor crftico 

de 0,485 


6.6 Criterio de falha de Bao-Wierzbicki 


O criterio de falha de Bao-Wierzbicki e fungao apenas do fator de 
triaxialidade, mas e representado por quatro fungoes diferentes de acordo com a 
faixa de triaxialidade. Este criterio de falha foi desenvolvido por Bao e Wierzbicki 
(2004) atraves da construgao do lugar geometrico de falha (fracture locus) pela 
realizagao e simulagao de varios testes experimentais incluindo testes em 
cisalhamento puro e testes uniaxiais de especimes em forma de cilindro com e sem 
entalhe, para o alurmnio 2024-T351. A fungao de deformagao de falha aparece 
esquematizada na segao 3.3.6 para o alurmnio 2024-T351, sendo que a mesma e 
ilustrada na Figura 55. 
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Figura 55. Deformagao de falha para o criterio de Bao-Wierzbicki 



Deformagao de engenharia (%) 

Figura 56. Resultado da simulagao com criterio de falha de Bao-Wierzbicki 

Apesar da formulagao deste criterio possuir uma meihor representagao da 
deformagao de falha para praticamente toda a faixa de triaxialidades, nao foi 
constatado um bom resultado numerico para a falha sob cisalhamento. Alem disso, 
uma grande oscilagao da curva de forga numerica foi verificada, conforme Figura 56. 

Apesar do resultado ruim para a previsao de falha por cisalhamento por parte 
deste criterio, em contrapartida pode-se dizer que o mesmo ja foi empregado com 
sucesso em simulagoes de impacto com falha, como em simulagoes de testes de 
Taylor com a verificagao de diferentes modos de falha do especime (Teng et al., 
2005). 
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Figura 57. Resultado da simulagao com criterio de falha de Bao-Wierzbicki 



Figura 58. Para o criterio de falha de Bao-Wierzbicki, duas trincas foram geradas na regiao do 

entalhe central 


6.7 Criterio de falha de Xue-Wierzbicki 


Este criterio considera o fator de triaxialidade e o angulo de Lode em sua 
formulagao, possuindo cinco parametros para a sua calibragao, os quais estao 
listados na Tabela 7. 


Tabela 7 - Parametros do criterio de falha de Xue-Wierzbicki para alumfnio 2024-T351 


c, 

c 2 

c 3 

c 4 

n 

0,87 

1,77 

0,21 

0,01 

0,153 


Na Tabela 7, n e o coeficiente de encruamento do material, extraido do ajuste 
da curva tensao-deformagao equivalente do material a fungao: 


v = K t e n p 


( 73 ) 
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onde crea tensao equivalente, K t e o modulo tangente e e p e a deformagao plastica 

equivalente. A superficie de deformagao de falha em fungao do fator de triaxialidade 
e angulo de Lode para este criterio e mostrada na Figura 59. 



Figura 59. Superficie de deformagao de falha para o criterio de falha de Xue-Wierzbicki. Extrafdo de 

Wierzbicki et al. (2005) 



Este criterio representou com uma razoavel acuracidade a falha no material, como 
pode ser verificado na Figura 60, alem de gerar uma falha simetrica em torno do 
entalhe 3. 
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Figura 61. Resultado da simulagao para o criterio de Xue-Wierzbicki 


6.8 Discussao dos resultados 


Com reiagao aos criterios baseados em um unico valor de calibragao, foi 
obtido um bom resultado para o criterio da maxima deformagao equivalente. Deve- 
se ressaltar que os resultados ruins dos criterios da maxima tensao cisalhante e do 
criterio de Cockcroft Latham dependem dos valores utilizados na calibragao, os 
quais geralmente sao calibrados para atender um maior numero de casos de testes 
experimentais de diferentes especimes. Com reiagao aos criterios de faiha baseados 
no fator de triaxialidade, obteve-se um bom resultado com o criterio de faiha de 
Johnson-Cook, apesar de sua simpiicidade. Para o criterio de faiha de Bao- 
Wierzbicki, somente a faiha em alta triaxialidade foi prevista com boa acuracidade. 
Assim, pode-se dizer que este criterio teve um resultado ruim levando em 
consideragao a sua complexidade. Para o criterio de faiha de Xue-Wierzbicki um 
bom resultado pode ser obtido. Deve-se destacar que o criterio de Xue-Wierzbicki e 
o unico criterio avaiiado aqui que tambem e fungao do angulo de Lode. 

Com reiagao ao formato de faiha dos especimes, esperava-se que ocorresse 
faiha simetrica. No entanto em alguns criterios de faiha verificou-se uma faiha nao 
simetrica em torno do entalhe central. Isso ocorreu devido ao fato de haver pequena 
rotagao em torno do entalhe central, o que gerava o imcio de duas trincas de lados 
opostos de forma nao simetrica, como iiustrado na Figura 54 para o criterio de 
Cockcroft-Latham. Para o criterio de Bao-Wierzbicki, inclusive, quase ocorreu a 
liberagao de um fragmento do entalhe central peio fato de haver ruptura nao- 
simetrica. 

Foi provado que o especime Bifailure e util para o benchmark de criterios de 
faiha. No entanto, uma limitagao do referido especime e o fato de ser possfvel 
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avaliar apenas duas triaxialidades diferentes. Desse modo, recomenda-se para 
trabalhos futuros a possibilidade do projeto de um especime como o Bifailure com 
falha em tres triaxialidades diferentes, provavelmente tornando diferentes as 
geometrias dos entalhes 1 e 2. 
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7 CONCLUSOES 


Apos um estudo detalhado da teoria de falha em materials ducteis, concluiu- 
se que a falha e influenciada por parametros como triaxialidade, nfvel de deformagao 
plastica e, de acordo com a literatura mais atual, angulo de Lode. Dessa forma, foi 
obtida a geometria de um corpo de prova especial de ensaio de tragao para 
benchmark de criterios de falha em dois estados de tensao distintos, o estado de 
tensao de alta triaxialidade e o estado de tensao proximo ao cisalhamento puro. 

O especime especial foi denominado Bifailure, sendo que a segunda 
geometria obtida para o mesmo resultou em baixa rotagao da regiao central do 
especime durante o teste, o que proporcionou a obtengao de estado de tensao 
proximo ao cisalhamento puro na regiao central, alem de medidas mais precisas 
com o extensometro. 

A fim de ser possfvel a comparagao de desempenho entre diferentes criterios 
de falha, uma subrotina especial foi desenvolvida para o programa comercial de 
elementos finitos LS-Dyna®. Verificou-se atraves da simulagao que criterios de falha 
mais simples, como o baseado na maxima deformagao plastica equivalente, 
puderam gerar resultados satisfatorios para o especime em questao. Inclusive, um 
criterio de falha tradicional como o de Johnson-Cook gerou resultados melhores do 
que criterios de fundamentagao mais complexa, como o de Bao-Wierzbicki e Xue- 
Wierzbicki. 

Acredita-se que os criterios de falha mais simples representam melhor 
situagoes de falha similares aquelas que ocorrem nos especimes utilizados para a 
calibragao. Por exemplo, o criterio de falha de Cockcroft-Latham e da maxima 
deformagao plastica equivalente possuem apenas um unico parametro de 
calibragao. Dessa forma, nao ha um valor crftico otimo comum a todos os tipos de 
falha. Se forem feitos varios testes, a diferentes condigoes, pode-se utilizar um valor 
que atenda razoavelmente bem a maior quantidade possfvel de situagoes. 

Para trabalhos futuros, propoe-se um reprojeto do especime Bifailure com 
tres entalhes de geometria diferentes, o que poderia proporcionar a falha em tres 
estados de triaxialidade diferentes ao inves de dois, como e o caso do especime 
atual. 
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APENDICE A - CALIBRAQAO DOS CRITERIOS DE FALHA DA 

mecAnica DO CONTINUO 
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Percebe-se que, em ensaios em corpos de prova, o fator de triaxialidade e 
angulo de Lode nao sao constantes durante todo o teste, conforme Figura 62. Dessa 
forma, aqui sera descrito um metodo de calibragao dos criterios de falha da 
mecanica do contfnuo que sao baseados na variavel cumulativa de dano: 


D =\ 



( 74 ) 



Figura 62. Obtengao do fator de triaxialidade medio para diferentes experimentos. Adaptado de Bao 

e Wierzbicki (2004). 

O processo de calibragao sera exemplificado aqui para o criterio de falha de 
Johnson-Cook, que apresenta a seguinte fungao de deformagao de falha e c 
(dependente apenas do fator de triaxialidade tj)\ 


£ c =D l +D 2 - exp {D 3 tj) 


( 75 ) 
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Aqui e necessario calibrar tres parametros: D,, D 2 e D 3 . Assim sao 

necessarios, no rnmirno, tres testes experimentais em diferentes valores de 
triaxialidade 77 . Neste caso, geralmente sao realizados testes em cisalhamento puro 
(77 proximo de zero) e testes de tragao uniaxial de especimes com e sem entalhe 
(altos valores de 77 ). Depois de realizados os testes experimentais, 0 proximo passo 
e realizar simulagoes em elementos finitos desses testes experimentais. Neste 
passo, e possivel empregar na simulagao modelos constitutivos que nao sao 
afetados pelo dano no material (sem perda de resistencia mecanica induzida por 
dano). Verifica-se na literatura que nesta fase sao empregados modelos de 
plasticidade que requerem uma curva tensao-deformagao verdadeira do material ou 
modelos de plasticidade como 0 de von Mises. O objetivo da simulagao numerica 
dos experimentos e extrair a deformagao plastica e os componentes de tensao ou 
deformagao da regiao onde ocorre ruptura dos corpos de prova, no momento de 
ruptura. Sabe-se, por exemplo, que para especimes com entalhe e possivel que a 
ruptura ocorra na regiao com 0 maior valor de triaxialidade, ou seja, no centra dos 
entalhes e nao no raio de curvatura dos mesmos. Assim, sabendo 0 local de 
ocorrencia de falha nos experimentos, nas simulagoes deve-se obter para 0 
elemento onde ocorreria a falha a deformagao plastica equivalente, e neste caso, 0 
fator de triaxialidade. No entanto, 0 fator de triaxialidade a ser obtido do elemento 
que falha e urn valor medio ao longo de todo 0 historico de carregamento do 
elemento durante a simulagao. Esse valor medio e dado pela formula: 



onde o m e a tensao hidrostatica e a e a tensao equivalente. Deste modo, para 0 

elemento onde ocorre a falha deve ser obtida a deformagao de falha e 0 fator de 
triaxialidade medio dado pela formula acima. A razao de ser utilizado urn fator de 
triaxialidade medio ao inves de urn fator de triaxialidade no instante de falha vem do 
fato de 0 equacionamento a seguir poder ser aplicado: 
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( 77 ) 


Para D = 1, temos: e f =s c {rj), o que possibilita a calibragao. Deste modo, 

obtidas as deformagoes de falha e as triaxialidades medias de ruptura para cada 
especime e possivel obter os parametros de falha utilizando a fungao de deformagao 
de falha £ c (jj) e o metodo dos mmirmos quadrados, e a partir daf montar um grafico 
que corresponderia ao lugar geometrico de falha, como o grafico da Figura 63. 



V m 


Figura 63. Obtengao do lugar geometrico de falha a partir dos valores de deformagao de falha e 

triaxialidade media obtidos numericamente. 


Um ponto importante a ser considerado e que para a calibragao e obtengao 
dos parametros de falha o fator de triaxialidade r/ na equagao de deformagao de 

falha e c torna-se o fator de triaxialidade medio rj m . 
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APENDICE B - SUBROTINA IMPLEMENTADA NO LS-DYNA® 




k 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

subroutine umat45 (cm,eps,sig,epsp,hsv, dtl, capa, etype, tt, 

1 temper,failel,crv) 
include 'nlqparm' 
include 'iounits.inc' 

dimension cm(*),eps(*),sig(*),hsv(*),crv(101,2,*) 
dimension aterms(4),rtr(3),rstava(7) 

real mult,bas,powp,effs,depi,A,B,n,qtrial,sigmay, phi, tol 
real ql,q2,q3,q4,q5,q6,dgama,g2,g,epbar,denom,ddgama,resnor 
real factor,R2G,R3G,davg,p,aj2,aux,auxl,aux2,aux3,epbar_i 
real eta,alfa 

integer nriter, mxiter, debug 
character*(*) etype 
logical failel 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 


Johnson-Cook solid visco-elastic-plastic material 
with isotropic hardening variables 
****** for BRICK elements ****** 


cm(1)=young's modulus 
cm(2)=poisson's ratio 
cm(3)=shear Modulus 

cm(4)=Hardening factor beta (Johnson-Cook parameter n) 
cm(5)=Tangent modulus (Johnson-Cook parameter B) 
cm(6)=yield stress (Johnson-Cook parameter A) 
cm(7)=Bulk Modulus 

cm(8)=C (Johnson-Cook material model parameter) 

cm(9)=eps0 (Johnson-Cook material model reference strain rate) 


cm(10)=Failure criteria 
(ifflag =0 
= 1 
=2 
=3 
= 4 
=5 


flag 

no failure, 

accumulated equivalent plastic strain 

maximum shear stress 

Johnson-Cook 

Wilkins 

Cockcroft-Latham 


=6 Bao-Wierzbicki 
=7 Xue-Wierzbicki 


★ 

cm 

(id 

= critical 

value for erosion (crit) 



★ 

cm 

(12) 

= constant 

for 

failure 

criteria 

(cfl) 




cm 

(13) 

= constant 

for 

failure 

criteria 

(cf 2) 



★ 

cm 

(14) 

= constant 

for 

failure 

criteria 

(cf 3) 



★ 

cm 

(15) 

= constant 

for 

failure 

criteria 

(cf 4) 



★ 

cm 

(16) 

= constant 

for 

failure 

criteria 

(cf5) - 

hardening exponent for 

★ 







Wierzbicki 

failure criterion 

k 

cm 

(17) 

=tolerance 

for 

Newton-! 

Raphson algorithm 

convergence 

k 

cm 

(18) 

=max numbe 

r of 

iterations for Newton-Raphson 

algorithm 


cm(19)=starting value for dgama, for Newton-Raphson algorithm 
cm(20)=debug variable, if equals 1 debug messages are printed 
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hsv(l)=total plastic strain 
hsv(2)=total pressure 
hsv(3)=von Mises stress 
hsv(4)=max shear stress 
hsv(5)=triaxiality (-hsv(2)/hsv(3)) 
hsv(6)=Lode parameter 
hsv(7)=cumulative damage 

hsv(8)=sum of (triaxiality)*(plastic strain increments) 
hsv(9)=sum of (Lode parameter)*(plastic strain increments) 

eps=strain increments to obtain elastic trial stress state 
ql-6=deviatoric stresses 
depi=incremental platic strain 


R0=0.E0 
Rl=l.E0 
R2=2,E0 
R3=3.E0 
R6=6.E0 
R1D2=R1/R2 
R1D3=R1/R3 
R2D3=R2/R3 
tol=cm(17) 
mxiter=cm(18) 
dgama=cm(19) 
debug=cm(20) 

★ 

* for failure 

k 

ifflag=cm(10) 
if(ifflag.ne.0.0) then 
crit=cm(11) 
cfl=cm(12) 
cf2=cm(13) 
cf3=cm(14) 
cf4=cm(15) 
cf5=cm(16) 
endif 

k 

* Set some material parameters 

* 


A=cm(6) 

B=cm(5) 
n=cm (4) 

C=cm(8) 
eps0=cm(9) 

★ 

* Shear modulus and auxiliary parameters 

k 

g2 =cm(l)/(Rl+cm(2)) 
g =RlD2*g2 
R2G=R2 *g 
R3G=R3*g 

k 

davg= (-eps (1) -eps (2) -eps (3) ) /3 . 

* Computing Hydrostatic Stress(Incremental) 

p=-davg*cm(1)/((R1-R2 *cm(2))) 



Updating Total Pressure 
hsv(2)=hsv(2)-p 

Computing Trial Stress 

sig(l)=sig(l)+p+g2*(eps(1)+davg) 
sig(2)=sig(2)+p+g2*(eps(2)+davg) 
sig(3)=sig(3)+p+g2*(eps(3)+davg) 
sig(4)=sig(4)+g*eps(4) 
sig(5)=sig(5)+g*eps(5) 
sig(6)=sig(6)+g*eps(6) 

Computing Deviatoric Stress 

ql=hsv(2)+sig(1) 
q2=hsv(2)+sig(2) 
q3=hsv(2)+sig(3) 
q4=sig(4) 
q5=sig(5) 
q6=sig(6) 

aj2=q4*q4+q5*q5+q6*q6+(ql*ql+q2*q2+q3*q3)/2. 


f(debug.eq.1) 

then 


write(59, *) 

'sig(1) 

= ', sig(1) 

write(59, *) 

'sig(2) 

= ', sig (2) 

write(59,*) 

'sig(3) 

= ', sig(3) 

write(59, *) 

'sig(4) 

= ', sig ( 4) 

write(59,*) 

'sig(5) 

= ', sig(5) 

write(59,*) 

' sig (6) 

= ', sig(6) 

write(59,*) 

' hsv (2) 

= ' , hsv (2) 

write(59,*) 

' a j2 = 

, a j2 


endif 


Calculation of equivalent stress and flow function 


qtrial=sqrt(R3*aj2) 
sigmay=A+B*(epsp**n) 
if(debug.eq.1) then 

write(59,*) 'initial_dgama = 


write(59,*) 
write(59,*) 
write(59,*) 
endif 


'R3G = ', R3G 
'initial qtrial = 
'initial sigmay = 


dgama 

qtrial 

sigmay 


Begin strain rate effect 

if(C.ne.RO) then 
dld=eps(1)+davg 
d2d=eps(2)+davg 
d3d=eps(3)+davg 
d4d=eps(4) 
d5d=eps(5) 
d6d=eps(6) 

ds=d4d* *R2+d5d**R2+d6d**R2 
Computing effective strain 

aux=R2 *(dld**R2+d2d**R2+d3d**R2+R2*ds)/R3 
if(aux.ne.RO) then 
effs=sqrt(aux) 
else 

effs=R0 
end if 

Computing strain rate 

if(tt.ne.RO) then 
effs=effs/dtl 
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Updating flow stress based on strain rate effect 
sigmay=sigmay*(Rl+C*log(effs/epsO)) 
end if 
end if 

Check for plastic admissibility 
phi=qtrial-sigmay 
if(phi/sigmay.gt.tol) then 

Plastic step: Apply return mapping - use Newton-Raphson algorithm 
to solve the return mapping equation 


epbar_i = epsp 

phi = qtrial - R3G*dgama - A - B*((epbar_i+dgama)**n) 
if(debug.eq.1) then 

write(59,*) 'Starting Newton-Raphson - phi=', phi 
write(59,*) 'initial epsp = ', epsp 
endif 

do 11 nriter=l,mxiter 

Compute residual derivative 
auxl=epbar_i+dgama 
aux2=n-Rl 
aux3=auxl* *aux2 
denom=-R3G-n*B*aux3 

Compute Newton-Raphson increment and update variable DGAMA 
ddgama=-phi/denom 
dgama=dgama+ddgama 
if(debug.eq.1) then 


'nriter = ', nriter 
'epbar_i+dgama = ', auxl 
'n-Rl = ', aux2 
'(epbar_i+dgama)**(n-Rl) = 


'denom = ' 
'ddgama = 

'dgama = ' 


denom 
. ddgama 
dgama 


write(59,*) 
write(59,*) 
write(59,*) 
write(59,*) 
write(59,*) 
write(59,*) 
write(59,*) 
endif 

Compute new residual 

epbar=epbar_i+dgama 
sigmay=A+B*(epbar**n) 
phi=qtrial-R3G*dgama-sigmay 
if(isnan(phi).eq.-l) then 

write(59,*) 'Not a number found!! phi = 
goto 997 
endif 

Check convergence 

resnor=abs(phi/sigmay) 
if(debug.eq.1) then 

write(59,*) 'isnan(phi) = ', isnan(phi) 
'epbar = ', epbar 
'sigmay = ', sigmay 
'phi = ', phi 
'resnor 


aux3 


phi 


resnor 


write(59,*) 
write(59,*) 
write(59,*) 
write(59,*) 
endif 

if(resnor.le.tol) then 
if(debug.eq.1) then 

write(59,*) '‘‘convergence achieved**' 

endif 

Update accumulated plastic strain 
epsp=epbar 

Update incremental plastic strain 
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depi=dgama 

Update stress components 

factor=Rl-R3G*dgama/qtrial 
sig(l)=factor*ql-hsv(2) 
sig(2)=factor*q2-hsv(2) 
sig(3)=factor*q3-hsv(2) 
sig(4)=factor*q4 
sig(5)=factor*q5 
sig(6)=factor*q6 
goto 997 
endif 

11 continue 

Reset failure flag and print warning message if the algorithm 
fails 

write (59,*) ' * CONVERGENCE FAILED**' 

else 

depi=R0 

endif 


* Updating accumulated plastic strain in history variable 

★ 

997 hsv(l)=epsp 

k 


post processor 


* von Mises stress 

vm=R2/R6*(sig(l)*sig(1)+sig(2)*sig(2) 
+ +sig(3)*sig(3)-(sig(1)*sig(2) 

+ +sig(2)*sig(3)+sig(3)*sig(1) ) ) 

+ +(sig(4)*sig(4)+sig(5)*sig(5) 

+ +sig(6)*sig(6)) 

vm=sqrt(R3*vm) 
hsv(3)=vm 


* Triaxiality 

if(vm.eq.RO) then 
hsv(5)=R0 
hsv(6)=R0 
goto 10 
else 

hsv(5)=-hsv(2)/vm 
end if 

★ 

* Calculation of principal stresses by diagonalization of stress tensor. 

* Finding the principal stresses for three-dimensional case means finding 

* the roots for the equation: s3 - AA*s^2 + BB*s - CC = 0 

k 

AA=sig(1)+sig(2)+sig(3) 

BB=sig(1)*sig(2)+sig(2)*sig(3)+sig(l)*sig(3)- 
+ sig(4)**2-sig(5)**2-sig(6)**2 

CC=sig(1)*sig(2)*sig(3)+R2*sig(4)*sig(5)*sig(6)- 
+ sig(1) *sig(5)**2-sig(2)*sig(6)**2-sig(3)*sig(4)**2 

aterms(1)=R1 
aterms(2)=-AA 
aterms(3)=BB 
aterms(4)=-CC 
call cubroot(aterms,rtr) 
if(rtr(3).gt.rtr(2)) then 
if(rtr(1).It.rtr(2)) then 
sigl=rtr(3) 



sig2=rtr(2) 
sig3=rtr(1) 

elseif(rtr(1) ,lt.rtr(3)) then 
sigl=rtr(3) 
sig2=rtr(1) 
sig3=rtr(2) 
else 

sigl=rtr(1) 
sig2=rtr(3) 
sig3=rtr(2) 
end if 
else 

if (rtr (1) .It.rtr (3)) then 
sigl=rtr(2) 
sig2=rtr(3) 
sig3=rtr(1) 

elseif(rtr(l).lt.rtr(2)) then 
sigl=rtr(2) 
sig2=rtr(1) 
sig3=rtr(3) 
else 

sigl=rtr(1) 
sig2=rtr(2) 
sig3=rtr(3) 
end if 
end if 

Ordering of principal deviatoric stresses 

rtr (1)=sigl+hsv(2) 
rtr(2)=sig2+hsv(2) 
rtr(3)=sig3+hsv(2) 

if(rtr(3).gt.rtr(2)) then 
if(rtr (1) .It.rtr (2)) then 
sl=rtr(3) 
s2=rtr(2) 
s3=rtr(1) 

elseif(rtr(1),lt.rtr(3)) then 
sl=rtr(3) 
s2=rtr (1) 
s3=rtr(2) 
else 

sl=rtr(1) 
s2=rtr(3) 
s3=rtr(2) 
end if 
else 

if (rtr (1) .It.rtr (3) ) then 
sl=rtr(2) 
s2=rtr(3) 
s3=rtr(1) 

elseif(rtr(1),lt.rtr(2)) then 
sl=rtr(2) 
s2=rtr(1) 
s3=rtr(3) 
else 

sl=rtr(1) 
s2=rtr(2) 
s3=rtr(3) 
end if 
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end if 

★ 

* Max shear stress 

hsv(4)=abs(sigl-sig3)/R2 

★ 

* Lode parameter 

hsv(6)=1.-(6./3.1416)*atan((2.*(s2-s3)/(sl-s3)-1.)/sqrt(3.)) 

k 

* Calculation of integrals for definition of average triaxiality and 

* average Lode parameter 

k 

hsv(8)=hsv(8)+hsv(5)*depi 

hsv(9)=hsv(9)+hsv(6)*depi 

k 

* Calculation of principal strain increments by diagonalization of 

* tensor made by components of strain increments. 

* Finding the principal stresses for three-dimensional case means finding 

* the roots for the equation: s3 - AA*s^2 + BB*s - CC = 0 

k 

AA=eps(1)+eps(2)+eps (3) 

BB=eps(1)*eps(2)+eps(2)*eps(3)+eps(1)*eps(3)- 
+ eps(4)**2-eps(5)**2-eps (6)**2 

CC=eps (1) *eps (2) *eps (3) +R2*eps (4) *eps (5) *eps (6) - 
+ eps (1) *eps (5) **2-eps (2) *eps (6) **2-eps (3) *eps (4) **2 

aterms(1)=R1 
aterms(2)=-AA 
aterms(3)=BB 
aterms(4)=-CC 
call cubroot(aterms,rtr) 
if(rtr(3).gt.rtr(2)) then 
if(rtr (1) .It.rtr(2)) then 
epsl=rtr(3) 
eps2=rtr(2) 
eps3=rtr(1) 

elseif(rtr(l) .It.rtr (3)) then 
epsl=rtr(3) 
eps2=rtr(1) 
eps3=rtr(2) 
else 

epsl=rtr(1) 
eps2=rtr(3) 
eps3=rtr(2) 
end if 
else 

if(rtr (1) .It.rtr(3)) then 
epsl=rtr(2) 
eps2=rtr(3) 
eps3=rtr(1) 

elseif(rtr(1),lt.rtr(2)) then 
epsl=rtr(2) 
eps2=rtr(1) 
eps3=rtr(3) 
else 

epsl=rtr(1) 
eps2=rtr(2) 
eps3=rtr(3) 
end if 
end if 


k k k 



k k k 


for different failure criteria 


if(ifflag.eq.0) return 
***Maximum equivalent plastic strain 
if(ifflag.eq.1) then 
hsv(7)=hsv(1) 

***Maximum shear stress 

elseif(ifflag.eq.2) then 
hsv(7)=abs(sigl-sig3)/R2 
***Johnson & Cook 

elseif(ifflag.eq.3) then 
eta=hsv(5) 

hsv(7)=hsv(7)+depi/(cfl+cf2*exp(cf3*eta)) 

***Wilkins 

elseif(ifflag.eq.4) then 

if((si.ne.RO).and.(s3.ne.RO)) then 
Al=s2/sl 
A2=s2/s3 

if(Al.gt.A2) then 
A=A1 
else 
A=A2 
end if 

hsv(7)=hsv(7)+depi*((R2-A)**cfl)/((Rl+cf2*hsv(2))**cf3) 
if((hsv(2).It.(-R1 /cf2)) .and.(hsv(1).gt.RO)) then 
hsv(7)=crit+Rl 
end if 
end if 

***Cockcroft-Latham 

elseif(ifflag.eq.5) then 
if(sigl.gt.RO) then 

hsv(7)=hsv(7)+depi*(sigl/vm) 
end if 

***Bao-Wierzbicki (for aluminum 2024-T351) 
elseif(ifflag.eq.6) then 
eta=hsv(5) 

if(eta.le.(-R1D3)) then 
hsv(7)=hsv(7)+depi/l.E10 
elseif((eta.gt.(-R1D3)).and.(eta.It.RO)) then 
hsv(7)=hsv(7)+depi/(0.1225*(eta+RlD3)**(-0.46)) 
elseif((eta.gt.RO) .and. (eta.It.0.4)) then 

hsv(7)=hsv(7)+depi/(1.9*eta**2.0-0.18*eta+0.21) 
else 

hsv(7)=hsv(7)+depi/(exp(-1.944*eta)) 
end if 

***Xue-Wierzbicki 

elseif(ifflag.eq.7) then 
eta=hsv(5) 

lode=(27./2.)*(sl*s2*s3)/(vm**3.) 

hsv(7)=hsv(7)+depi/(cfl*exp(-cf2*eta)-(cfl*exp(-cf2*eta) 

+ -cf3*exp(-cf4*eta))*((1.-lode**(1./cf5))**cf5)) 
end if 

if (hsv(7).ge.crit) then 
sig(1)=0.0 
sig(2)=0.0 
sig(3)=0.0 
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sig(4)=0.0 
sig(5)=0.0 
sig(6)=0.0 
failel=.true. 
end if 
return 
end 


10 



